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Abstract: We study decoupling in FRW spacetimes, emphasizing a Lagrangian descrip- 
tion throughout. To account for the vacuum choice ambiguity in cosmological settings, 
we introduce an arbitrary boundary action representing the initial conditions. RG flow 
in these spacetimes naturally affects the boundary interactions. As a consequence the 
boundary conditions are sensitive to high-energy physics through irrelevant terms in the 
boundary action. Using scalar field theory as an example, we derive the leading dimension 
four irrelevant boundary operators. We discuss how the known vacuum choices, e.g. the 
Bunch-Davies vacuum, appear in the Lagrangian description and square with decoupling. 
For all choices of boundary conditions encoded by relevant boundary operators, of which the 
known ones are a subset, backreaction is under control. All, moreover, will generically feel 
the influence of high-energy physics through irrelevant (dimension four) boundary correc- 
tions. Having established a coherent effective field theory framework including the vacuum 
choice ambiguity, we derive an explicit expression for the power spectrum of inflationary 
density perturbations including the leading high energy corrections. In accordance with 
the dimensionality of the leading irrelevant operators, the effect of high energy physics is 
linearly proportional to the Hubble radius H and the scale of new physics I = 1/M. 
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1. Introduction 



String theory provides a fundamental framework to describe physics at the highest energy 
scales. Yet, the details of transplanckian physics have completely eluded us so far. For- 
tunately, the notion of decoupling allows us to understand low energy phenomena despite 
our ignorance of physics at very high energies. Renormalization Group (RG) flow teaches 
us that the effects of high energy physics can be captured by only a finite number of rel- 
evant couplings in the low energy theory. In flat spacetime, the decoupling between high 
and low energy physics is well established. However, for quantum field theories in curved 
space and in FRW universes in particular, decoupling is not so clearcut. In cosmological 
spacetimes high energy scales are redshifted to low energy scales via cosmic expansion. 
This connects high and low energy physics through unitary time evolution in addition to 
the dynamics. Decoupling, specifically in the inflationary context, is of great importance 
to upcoming cosmological precision experiments. All current physical scales would origi- 
nate from transplanckian scales at the onset of inflation, if inflation lasted longer than the 
minimal number of e-folds. Conceivably, then, signatures of Planck scale physics (stringy 
or other) could show up in cosmological measurements |l|, ||, ||, ||, ||, ^ [7] . This possibil- 
ity whether glimpses of transplanckian physics can be observed in the cosmic microwave 
background (CMB) radiation || is determined by the strength with which transplanckian 
physics decouples. Remarkably, such effects are potentially observable, but only if the 
transplanckian physics selects a non-standard initial state H]. 1 Other high energy ef- 
fects are generically too small Q (with the exception of the higher dimensional operators 
identified in |||). More recently, explicit examples were presented to illustrate that the 
integrating out of a massive field could result in a non-trivial initial state, offering both 
a proof of principle that transplanckian physics may be observable, and suggesting that 
decoupling is more subtle in expanding universes Q. 

In this article we would like to clarify the connections between vacuum/initial state 
selection and decoupling in a fixed FRW background (we ignore gravitational dynamics 
throughout). In cosmological settings, i.e. in a spatially homogeneous and isotropic uni- 
verse, the size of the scale factor yields a preferred time coordinate, and as a consequence 
a Hamiltonian approach has become standard ||. In contrast to the Hamiltonian point 
of view which emphasizes the dynamical evolution, a Lagrangian point of view emphasizes 
the symmetries and scaling behaviour relevant to physical processes (see e.g. @, 0, [H]]). 
It is therefore the natural framework for a Wilsonian RG understanding of decoupling of 
energy scales and relevant degrees of freedom determined by symmetries. 2 However, a La- 
grangian or an action by itself is insufficient to determine the full kinematic and dynamic 
behaviour of quantum fields. One must in addition specify the boundary conditions. This 
corresponds to the choice of initial or vacuum state in the Hamiltonian language. The 

1 The nomenclature 'non-standard vacuum' state is also used. Strictly speaking there is no clear vacuum 
state in an FRW universe. In an abuse of language, we use vacuum and initial state interchangeably. 

2 Wilsonian RG in effect explains why (non-gravitational) physics works. Its success strongly suggests 
that the same principles are at work in quantum gravity and that general relativity is the low energy 
effective action relevant at scales below Mpianck (for a nice review on general relativity as an effective field 
theory see fff). String theory, in particular, is an explicit manifestation of this idea. 
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question directly relevant to the window on transplanckian physics provided by inflation 
is therefore which boundary conditions to impose on the fields. To preserve the symme- 
tries of the Lagrangian a subset of all possible boundary conditions is often only allowed. 
With enough symmetry, e.g. Minkowski QFT, the choice may in fact be unique. FRW 
spacetimes have less symmetry and it is a priori not clear, what the natural or correct 
boundary conditions are. What we will explain in section ^ is that no matter which choice 
of boundary conditions is made in the full quantum theory, RG-flow in the effective low 
energy action will generically change these conditions. In particular high-energy physics 
will affect the boundary conditions through irrelevant corrections, which we derive. We 
apply these results in section [| to the computation of the power spectrum of inflationary 
density perturbations. The leading irrelevant correction to the boundary conditions is of 
dimension four, and we therefore find that the power spectrum is subject to corrections 
of order H/M with M the scale of new physics. This is in accordance with earlier pre- 
dictions that transplanckian effects are potentially observable [||, 0. Importantly, we are 
able to derive this result purely within the framework of Wilsonian effective field theory. 
This makes our answer predictive both in the sense that the parametric dependence of 
inflationary physics on high-energy is now manifest, and that the strength is computable 
in any theory where the high energy physics is explicitly known. Because our results are 
derived within the context of effective field theory, they provide a settlement to the de- 
bate [||, [| ||, 12] whether H/M corrections are consistent with decoupling arguments. We 



conclude with an outlook where we will briefly comment on the relation of our results to 
consistency issues regarding (non-trivial) de Sitter invariant vacua known as a-states. We 
will, however, begin with a summary, lest the trees obscure the forest. 

1.1 Summary of our results 

Any boundary conditions one wishes to impose can be encoded in a boundary action. This 



is even true for the Minkowski vacuum (section |2.3p . It has long been known that the 
couplings in such a boundary action are renormalized at the quantum level. Equivalently, 
a Wilsonian approach to the effective action ought to result not only in a renormalization 
of the boundary couplings, but also in the generation of irrelevant boundary operators. 
Consider, for example, a two scalar field model with a mass separation M x S> and 
boundary and bulk interactions S znt = — J gx4> ~ §ix4>- This is exactly solvable, and 
upon integrating out x> permitted when the cut-off scale A <C M x , one generates the 
boundary interactions 

A, A, 

We will describe and review the Wilsonian effective action for theories with a boundary, 
including this example, in section ^. 

The issue of (boundary) Wilsonian decoupling is relevant to our understanding of 
cosmology. In an expanding universe, there is no unique vacuum state. In the Lagrangian 
language, this translates to a lack of knowledge of the appropriate boundary conditions. 
Recall that any boundary conditions, including the 'Minkowski' ones, can be encoded in 
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a boundary action. Wishing to emphasize the Lagrangian viewpoint, where the study of 
decoupling is most natural, we add a boundary action with free parameters at a fixed but 
arbitrary time to. 

Our limited understanding of high-energy physics in the very early universe can thus be 
accounted for by the inclusion of a boundary action in a cosmological effective Lagrangian. 
Whichever boundary conditions we choose this boundary action to encode, they will be 
subject to renormalization. In particular, the details of the high-energy physics, which has 
been integrated out, will be encoded in irrelevant corrections to the boundary action. For 
Z2 symmetric scalar field theory the leading irrelevant boundary operators (that respect 
the homogeneity and isotropy of FRW cosmologies) are 



qirr.op. 
bound 



drx 



- A 4 



(1.2) 



where d n is the normal derivative. These operators are of dimension four — one di- 
mension higher than the boundary measure — and describe corrections of order \k\/M 
plus a boundary four-point interaction. For the momentum range of interest to the CMB, 
\k\ ~ H, where H is the Hubble parameter, the quadratic operators scale as H/M and they 
are therefore the primary candidates for witnessing consequences of high-energy physics in 
cosmological data. The leading bulk operator is of order H 2 /M 2 and is generically beyond 
observational reach Q. Computing the inflationary perturbation spectrum in a de Sitter 
background, including the corrections to Bunch-Davies boundary conditions due to the 
irrelevant operators (|1.2j), we find 



pdS 

BD+ 
irr.op. 



with 



7T 
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(1.3) 



KBD 



d-l + 2v 



H 



\k\ H v+ i(y 
a H u (y ) 



(1.4) 



where H u (yo) are Hankel functions at yo = \k\/a(r]o)H whose index v(m 2 ) depends on the 



mass m . Crucial in our exposition will be the proof (section 2.2) that, despite appearances 



this expression does not depend on the location of the boundary action yo- Only the 
meaning of the initial conditions matters, not where they are imposed. 

Eq. ( |1.3| ) is our main result. Having translated the cosmological vacuum choice am- 
biguity into an arbitrary boundary action, we conclude based on Wilsonian decoupling 
that the leading irrelevant operators in FRW field theory are boundary operators at order 
H/M. Using optimistic but not untypical estimates of H ~ 10 14 GeV and M ~ 10 16 
GeV (string scale), new (transplanckian) physics will generically affect the standard pre- 
dictions of inflationary cosmology at the one-percent level. Conversely, CMB observations 
with an accuracy of one percent or better can potentially measure effects of transplanckian 
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physics. Only for very special choices of initial conditions and transplanckian physics will 
this correction be absent. 

We further identify the boundary conditions corresponding to several cosmological vac- 
uum choices including the generalization of the "Minkowski-space" boundary conditions 
(sections |2.3| and |3.1| ). In the Wilsonian effective Lagrangian description it is clear that 
no vacuum is preselected by a consistency condition. Any boundary condition encoded by 
relevant operators is consistent, in the sense that the Minkowski space stress tensor coun- 
terterm generated with the appropriate boundary conditions will render the cosmological 
stress tensor finite as well (section ||). Backreaction is always under control. Which cos- 
mological boundary conditions are the right ones to impose, requires just physical input, 
as it should be. 



2. Decoupling in theories with a boundary: a review 

The study of field theories is primarily concerned with Minkowski backgrounds, with the 
symmetry-compatible boundary conditions that the fields vanish at infinity. 3 Actions which 
contain explicit boundary interactions, however, have been studied in the past [13, 14, 15 



|lq| , and are receiving renewed attention (see e.g. |17|, [L8|, [L9|, |20, 21, ^2[). One can use such 
boundary interactions to enforce whichever boundary conditions one wishes. Consider, for 
example, scalar A</> 4 theory on a semi-infinite space 4 

S bu ik = / d 3 xdy - -{d^f - —tf - , (2.1) 

Jyo<y<oo ^ Z 4. 

with the following boundary interactions added 

Sboundary = j> d 3 X~ ^(pd n (p - ^(f) 2 . (2.2) 

Here d n = d y is the derivative normal to the boundary. Expanding the action to first order 
in <j) + 5(j), we find the usual equation of motion 

SS bulk = J d 3 xdy 8<j> (n<j> - m 2 4> - , (2.3) 



plus the boundary conditions 



<W,„, / ~ / <f':r - <V> ( ^cW> + kA - | ft^ . (2.4) 



3 One may alternatively think of Minkowski space field theory as defined on a (infinite volume) torus 
("putting it in a box"), which has no boundary at all. 

4 We choose Lorentzian + + +— signature throughout the paper. Working with effective actions, we 
implicitly assume that all results can be obtained by a Wick rotation from Euclidean space. Depending 
on whether the boundary under consideration is spacelike or timelike relative signs and factors of i will 
appear. Our focus will be on spacelike boundaries in particular since those have a natural interpretation 
as initial states in a Hamiltonian description. We discuss details relating to the signature of time and the 
Wick rotation of timelike to spacelike boundaries in appendix [b| 
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If we insist that the variations S<p are arbitrary and do not vanish on the boundary (which 
would correspond to imposing Dirichlet boundary conditions), it appears that [A must 
vanish for consistency. As we will see shortly, however, renormalization can produce coun- 
terterms proportional to fj, and a more correct point of view is that <j) can be discontinuously 
redefined on the boundary pq] , together with a redefinition of the couplings which absorbs 

4>(x, y) <j)(x, y) + a6(y - y)4>{x, y ) , 

/ ( a2 \ f a<2 /ua 2 \ 2/i 

K = k + k I a + — I +<5(0) I — - fxa - — I , a = _ ■ (2.5) 

This field redefinition can be interpreted as a shift of the boundary value of <j) to the correct 
saddlepoint. 6 That this is the correct interpretation follows from the fact that we can also 
treat fi perturbatively as an interaction. A Feynman diagram computation will then yield 
an effective action with coupling k' 7 After this 'renormalization' the boundary term from 
partial integration is canonical 

SS bound = j> d 3 x - 8(j)d n (j) - k54>4> (2.6) 

which vanishes when 

d n ct> = -K'<f> . (2.7) 

We see that the (renormalized) value of k determines the boundary condition. For k = 
we have Neumann boundary conditions, for k = ±oo the (particular) Dirichlet boundary 
condition 4>(x, yo) = 0, and for finite k a mixture of the two. All possible (linear) boundary 
conditions are recovered. This is comforting as there are no other terms of order 4> 2 
compatible with the symmetries. In fact, the boundary action Sbound is the most general 
one we can write down, if we limit our attention to relevant operators 8 and require (for the 
sake of simplicity) that the action is also invariant under the bulk Z2 symmetry <p «-> —(f). Of 
course, for a second order PDE one needs two boundary conditions. The other comes from 
the second boundary of integration. In the example above this is y = 00. See appendix [A| 
for details. 

J Here 6(y) is the step function, with 8(0) — 1/2 and d y 9(y) — S(y). Recall that this distribution is of 
measure zero, i.e. dy8(yo — y)f(y) = 0. Of the bulk terms only the kinetic term is therefore affected by 
the shift. Also note that 5(y - yo)f(y) = 5/(2/0)- 

One can also find a redefinition of the type <j>'(y) = <j>{y) + ct6(yo — y)<f>(y), which is the correct one from 
the point of view of coarse graining and the distributional definitions for 9(y) and S(y). Interestingly, the 
redefinitions required are the same. 

6 When counterterms of the form <j>d n ^> are required for renormalization, this shift of the background 
value for <j> is thus a boundary analogue of the Coleman- Weinberg phenomenon. 

7 A perturbative comparison with Feynman diagrams, which we perform in appendix [d|, also explains the 
delta function at zero argument. It serves to make all distributions conform to the bare boundary condition 
d n 4> = — kc/). 

8 We assume that the initial state encoded by the boundary action Sbound has no intrinsic size, i.e. 
a dimensionful scale. We are ultimately interested in vacuum-like initial conditions in cosmology. This 
restriction to scale-less initial states is therefore a natural one. 
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RG arguments then tell us, that in a bounded space the terms in the boundary action, 
even if they were not present at the outset, would be generated as counterterms. They are 
necessary for the consistency of the theory. Let us show this explicitly. Suppose we start 
with Neumann boundary conditions: k initially vanishes. By the method of images, the 
Neumann propagator equals 9 

f d?k dk e ikx ( Xl ~ X2 ^ (e iky ( yi ~ y2 ^ + e iky (~ yi ~ y2+2y °}) 
G N ( Xl ,y 1]X2 ,y 2 ) = -ij -^JL k 2 + k 2 + m 2 ■ ( 2 - 8 ) 

We will choose to regulate our theory by multiplying the propagator by a regulating func- 
tion ^"(D/A 2 ) = exp(— k 2 /A 2 ) [24]. This makes the path integral well defined and cleanly 



separates out the ultraviolet divergences. The one-loop seagull graph then evaluates to 




{<fi(xi,yi)4>(x 2 ,y2))i-ioop 



-iX 



G N (x 1 ,yi;x 1 ,yi)5 3 (x 1 - x 2 )5(y 1 - y 2 ) 



~ A ^;l,2^;l,2 

4(2^) 4 



d*k 



e a 



k 2 + m 2 



+ / d 6 k x dk 



Jk y {~2y+2y )-K I 

" k 2 + k 2 + m 2 



(2.9) 



The first term is the usual bulk \<fi A divergence of the two-point function. The second 
term, however, is a newly divergent term, and quite obviously a direct consequence of the 
boundary conditions. Evaluating this term in more detail, we find 



l—loop 



XS l,l,2 S y; 



1,2 



, vr 5 / 2 AeT 
4(2tt) 4 ' 

XA5 3 J( yi - y 2 )6( yi 



A 

A Jo 

y °)|A 2 >m 2 



(2.10) 



Note that the new divergence is entirely located on the boundary. The last step utilizes one 
of the more common distributional definitions of the Dirac-delta function (before doing the 
finite integral over s). Recalling the coarse-graining steps underlying RG-flow, it should 
come as no surprise that the delta- function localization appears in a distributional limit. 
This simply reflects that our spatial resolution decreases under RG-flow, and the precise 
location of the boundary becomes fuzzy. 

That the divergence is concentrated solely on the boundary (in this distributional 
sense) is reassuring. Bulk UV-physics should be unaffected by the presence of a boundary. 
It is precisely the breaking of Lorentz invariance due to the presence of the boundary that 
is responsible for the new divergence. By necessity it must then appear in the same sector 
of the theory that was responsible for the symmetry-violation in the first place. 

9 Our domain of interest y £ [yo, oo) is semi-infinite. Hence k v is a continuous variable. 
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To make the theory finite, we therefore need to add a boundary counterterm of the 
type 10 

S£SS = f <?x £ 2 (m 2 /A 2 ) (^) 2 ■ (2.11) 

with £(m 2 /A 2 ) chosen such that it cancels the divergence in eq. (2.10). This result is of 
course expected (in part) purely on dimensional grounds. 

The necessity of this counterterm has serious implications, however. Recalling the 
results from the first half of this section, we see that the boundary conditions change under 
RG-flow. In order to reproduce the same physics in a theory with a different cut-off, we 
not only need to change the vertices, but also the boundary conditions. (More precisely, 
to maintain a given physical renormalized boundary condition K ren we need to change the 
bare coupling k.) Of course, this counterterm is scheme-dependent. The beta-functions 
at one loop on the other hand are scheme- independent, and we can extract the generic 
behaviour of the boundary conditions from them. We find that as we change the scale, the 
boundary conditions change under RG-flow as 



~ <9A 



C A^W + 0(A 2 ) • (2.12) 

m 2 / A 2 fixed 77 ' 



with £ 2 > 0. This may seem surprising, but it does not go against the lore that boundary 
conditions are determined by physical conditions, and not by dynamics. It is worthwhile 
to repeat that what the RG-scaling of the boundary conditions says, is that in a cut-off 
theory, under a change of the cut-off, one reproduces the same physics when one changes 
the boundary conditions according to eq. ( |2.12 ). 



2.1 Boundary RG fixed points and 'vacua' 

A natural question to ask is what the endpoints of boundary RG-flow are. The explicit 
dimensionality of the coupling k already betrays the answer. In the deep IR, when \p\ <C A 
(A — » oo effectively; m = //A), k blows up, and the boundary conditions tend to the 
special Dirichlet boundary condition <fr(x,yo) = 0. Physically this is easily understood in 
Wilsonian RG language. The moment the cut-off restricts the momentum scales \p\ to be 
smaller than m (A ~ m), all modes freeze out and the theory ceases to be dynamical. 
Hence the field <f> 'vanishes', and must be Dirichlet. 

Dirichlet conditions thus form a trivial fixed point of RG-flow. This is easily visible. 
When (j) strictly vanishes on the boundary, simply no counterterms are possible. Both 
terms <f> <pd n (j) and <f> 4> 2 vanish. For completeness, were one to repeat the computation 
eq. ( [2.S| ) for Dirichlet conditions, the difference is that the propagator now has a relative 
minus sign. As a consequence, the bulk divergence cancels the boundary divergence at 
y = uq. Eq. ( |2.9| ) shows this clearly. In effective field theory the distinction between the 
fuzzy boundary and the bulk disappears in the deep IR limit, which explains why we can 
no longer treat bulk and boundary singularities separately when the boundary conditions 
become Dirichlet. 



3 Since the 'bare' boundary conditions are Neumann, this is the only type we can add. 
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When the boundary is spacelike and represents initial conditions in time, the induced 
changes in the boundary conditions due to RG-flow have a natural description in the 
Hamiltonian language of states. Under coarse graining the original state gets screened by 
vacuum polarization. In the low-energy effective theory, the correct state to use is a dressed 
version of the original state. If we take this picture further, we can deduce the boundary 
conditions which correspond to the vacuum. If the vacuum is the 'empty' state, then it 
ought not to become dressed under coarse graining. Translating back to the Lagrangian 
language, this means that the corresponding boundary conditions will not suffer from 
renormalization. Hence a vacuum in the Hamiltonian language should correspond to a 
fixed point of boundary RG-flow. 11 

2.2 Freedom of choice for the boundary location 

What will be of fundamental importance to us, is that the location of the boundary is 
arbitrary. The introduction of a boundary action at yo is a way to encode the initial 
conditions at the level of the action, but it does not necessarily mean that there is a 
physical object or obstruction at y = yo. It is simply a translation of the statement that a 
second order PDE needs two boundary conditions, but at what location one imposes those 
conditions is irrelevant. Of course, if one imposes the boundary conditions at a different 
location, they will not in general be of the same form as the original initial conditions. If 
one changes the location y one must change the value of k to keep the physics unchanged. 
A symmetry is therefore present between the location yo and ^- 12 To show this explicitly, 
choose a basis <p + (k,y), ip-(k,y) = ip* + {k,y) for the two independent solutions of the 
kinetic operator. In terms of this basis, the linear combination which obeys the boundary 
condition d n ip(yo) = —Kip(yo) is 

<p hK (k,y) = <p + (k,y) + b K (k)<p-(k,y), b K {k) = - ^ + '° ± ^+'° , (2.13) 

Here the subscript means that the quantity is evaluated at the boundary yo- Obviously 
if b K stays the same, physics stays the same. This allows us to derive a symmetry relation 
between the value k and the location yo- Under a constant shift of the boundary 5ip = 
£,d n (p = £d y ip and a simultaneous change 5k, b K changes as 13 



5b K = -Z 
—8k 



Kd n ip +fi + d%ip +j0 Kip +fi + d n tp +fi 2 
r~« 7 ^ («»n¥>-,o + d n ip- ) 



tup- ft + d n ip- fi (Kip- fi + d n ip- fi )< 



i<P-,o) 



(2.14) 



11 Presumably this is a UV-fixed point. Exciting the vacuum to a state, i.e. deforming away from the 
fixed point, reinstates RG-flow. The excitation, however, should not disappear in the deep IR. Hence the 
dressing of the state due to coarse graining leads one away from the vacuum. Of course to study boundary 
RG-flow, one needs an interacting theory. Any state in a free theory is a trivial fixed point of boundary 
RG-flow. 

12 This is not a true symmetry of the action. Because the coupling constant k changes, it is an isomorphism 
between families of theories. This is analogous to general coordinate invariance of the target space manifold 
in non-linear sigma models. 

13 Note that b K depends on the basis choice tp±, but k does not. 
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Demanding that 5b K vanishes, one finds the change in n necessary to keep physics unchanged 
under a change of the location of the boundary. This shows explicitly that this location is 
arbitrary. 

2.3 Minkowski space boundary conditions 

Minkowski space formally does not have a boundary of course. The arbitrariness of the 
location of the boundary, however, suggests that we should be able to treat it in a similar 
way. This is not quite manifest because, to stay within the framework of effective field 
theory, k must remain an analytic dimension one operator in the spatial momenta. The 
symmetry ( |2.14 ) is subject to this condition. The harmonic oscillator boundary conditions, 



constructed here to yield physics equivalent to unbounded Minkowski space physics, will 
be consistent with this requirement. To find these conditions suppose the boundary is 
a fixed time slice. We can then take a cue from the Hamiltonian formalism. Minkowski 
boundary conditions should correspond to choosing the standard Minkowski vacuum in the 
Hamiltonian picture. By definition this is the state annihilated by the lowering operator 
of each spatial momentum mode k x (in the free theory). 



a^|0)=0 (n^-iuj(k,m)^A |0) = , w(jfe, m) = \J 'k 2 + m 2 . (2.15) 

The canonical momentum conjugate to irk = do<pk is precisely the normal derivative to the 
fixed time slice. This suggests that we should choose the spatial momentum dependent 
boundary conditions pf] 



d n cj)\y = y = iy k 2 + m 2 (j>\y= yo — ► k = -iy k 2 + m 2 . (2.16) 

This boundary condition descends from the 'higher derivative' operator <f (j)^J d 2 — m 2 (f). 
But, as k has canonical dimension one, there is no new scale associated with this higher 
derivative term. Note that k is purely imaginary. This is a consequence of imposing the 
boundary condition at a fixed time. Wick rotating from a spatial boundary with real n 
generates a factor of i in the boundary condition d(p = —n<p. We provide details behind 
this naive argument in appendix [B|. We show there that all correlation functions will be 
analytic in the boundary coupling k, as is usual in effective field theory. We are therefore 
instructed to treat k as real throughout all steps of the calculation, and only substitute its 
imaginary value at the end. 

This momentum dependent choice of boundary conditions indeed ensures that the 
theory reproduces Minkowski space dynamics. For an arbitrary k the Green's function is 
(see eq. ( 2.1 3| ) , and recall that y parametrizes a timelike direction) 



. . . f d^kdky^ l e ^ ik y -K C . 

G.iw,*,*) = -,j -^-f g-qw-fc ■ (2 ' 17) 

where we have included the ie term. The second term, at first sight, negates equivalence 
with the Minkowski propagator 

f d 3 kdk e^i-^H^Ms/i-J^) 

GrMink = — i f ~, ~ > (2-18) 

J (2?r) 4 k 2 - k 2 + m 2 - ie 
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The coefficient k, however, is precisely chosen such that on shell the second term vanishes. 14 
By unitarity, the theory with n = —icu(k, m) is then the same as the Minkowski space 
theory. We can see this explicitly by performing the integral over k y . Doing so returns the 
standard Minkowski propagator in Hamiltonian form 

/^3^, e ik(xi—X2)—iu{k,m)(yi—y2) 
2^ 0(y 1 - y 2 ) + (y 2 ^ yi) , (2.19) 

which shows that the second term really is spurious. Indeed, this choice of k removes the 
pole in the second term, which means its contribution to any physical quantity disappears. 

We still have an official boundary at yo of course, even though the specific boundary 
conditions ( p. 16 ) ensure that it has no effect on physical amplitudes. The situation de- 



scribed here, is familiar from electrodynamics. 15 We have chosen an interface at yo where 
the dielectric properties happen to be the same for both materials. The transmission coef- 
ficient is therefore 100% and the wavefunction behaves as if the interface is not there, i.e. 
the interface is completely transparent. 

2.3.1 Minkowski boundary conditions and RG-flow 

Classical physics is indeed insensitive to a completely transparent interface. Is the quantum 
physics as well? In other words does the fact that the off-shell propagators appear to differ 
become relevant at the loop level? The answer is obviously no in perturbation theory. The 
cancellation of the pole by the specific 'Minkowski' choice for k means that in any integral 
the contribution of the second term vanishes. Hence the Minkowski boundary conditions 
do not get renormalized. They are a fixed point of boundary RG-flow exactly as befits the 
boundary conditions corresponding to a true vacuum. The reason why this is so is clear. 
The choice KMink = —ioj{k,m) is precisely the one that restores the Lorentz symmetry 
naively broken by the introduction of a boundary. Counterterms are forbidden to appear 
for they would break the reinstated Lorentz symmetry. 

2.4 Wilsonian RG-flow and irrelevant operators 

Quite generically therefore the boundary conditions of a quantum field theory are affected 
by RG flow, unless they are protected by a symmetry. Integrating out high energy degrees 
of freedom necessitates a change in boundary conditions to reproduce the same physics in a 
low-energy effective description of the theory. Decoupling then ensures that the low-energy 
theory remains predictive: the effects of high-energy physics are primarily encoded in a 
small set of relevant operators with universal scaling behaviour independent of the details 
of the high-energy theory. Subleading corrections of an energy expansion are by definition 
captured by irrelevant operators. These encode the specifics of the high-energy completion 
of the theory. 

One of our best hopes to detect the properties of high energy physics beyond the Planck 
scale is in a cosmological setting. The tremendous cosmological redshift during inflation 



14 The second term only vanishes for the domain 8(yi + yi — 2yo). Since our domain of interest is y > yo, 
this is always true. 

15 Except that this boundary is spacelike, which is why we can in fact relate it to a choice of initial state. 
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may bring the consequences of such irrelevant operators within reach of experimental mea- 
surements. This exciting opportunity has been a preeminent question in recent literature. 
In section || we shall show that the irrelevant boundary operators discussed in this subsec- 
tion are responsible for the leading effects of high-energy physics in cosmology, appearing 
generically at order H/Mpi anc k. The leading irrelevant operators for the bulk theory have 
long been known and their consequences for cosmological measurements are discussed in 
Q. However, it is well known that quantum field theory in cosmological settings suffers 
from a vacuum choice ambiguity. In the Lagrangian language this corresponds to a choice 
of boundary conditions. As we have just seen, we can parametrize this ambiguity in the 
cosmological vacuum choice by adding an arbitrary boundary action <f> ncf) 2 . Whichever the 
value of k may be, the influence of high-energy physics will be encoded in the irrelevant cor- 
rections to the boundary action. For that reason, we devote this section to a determination 
of the leading irrelevant operators on the boundary. Earlier studies have indeed indicated 
it is only (irrelevant) changes in the boundary condition which can have observable effects 
in measurements. Due to the symmetry constraints on the action the consequences of bulk 
irrelevant operators are just too small to be detectable. Our aim here is to provide a solid 
foundation for these earlier results. 

One can make a straightforward guess as to what the leading boundary irrelevant op- 
erators are, insisting on locality, compatibility with the Z2 symmetry, and SO (3) rotational 
invariance on the boundary. 16 They are the dimension four operators: 



d 6 x cf , * d A x d l 4>di(j) , <b d 6 x d n (j)d n (j) , <b d 6 x 4>d n d n <t> . (2.20) 
y=yo Jy=yo Jy=yo Jy=yo 

Note that the breaking of Lorentz invariance on the boundary distinguishes normal and 
tangential derivatives, and that normal derivatives cannot be integrated by parts. Varying 
<j> infinitesimally, the latter two will generate normal derivatives on the variation d n 6<p. To 
restore the applicability of the calculus of variations, one needs to perform a discontinuous 
field redefinition and adjustment of the couplings similar to fl2.5|) . (We do so in appendix 
g) In this sense, all physics can be captured by the first two irrelevant operators. However, 
for tractability we will treat all four operators perturbatively and on the same footing. We 
will see in section ||] that these operators will lead to corrections of order H/Mpi anc k to 
inflationary density perturbations, as predicted by the studies P]. Here we will give an 
explicit example where high-energy physics induces two of these dimension four irrelevant 
boundary operators. 

Tree-level diagrams exchanging a heavy field are the natural candidates for producing 
higher derivative corrections under RG-flow. We therefore add a scalar x to the theory 
with mass M x > A, to represent the high energy sector whose influence we will deduce. 
The only communication between the field \ an d <t> will be through the 'flavor-mixing' bulk 
and boundary couplings 

Sugh = ~ f d 3 xdy g X 4> ~ j> d 3 x 1X 4> , (2-21) 



16 These symmetry constraints follow from the assumption that the initial state has no intrinsic dimen- 
sionful parameter. See footnote 
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and x w i n have no other bulk or boundary (self)-interactions. Because the mass of x is 
higher than the cut-off, it will not appear as a final state, and in this simple model we 
can integrate it out explicitly. Its influence on the low-energy effective A</> 4 theory is only 
through tree-level mass oscillation graphs and a boundary reflection. Treating the couplings 
g and 7 as perturbations — hence the propagator for x wm have Neumann boundary con- 
ditions — consider the tree level correction to (4><f>) represented by the following Feynman 
diagram and its effective replacement. 



(2.22) 



Here wiggled lines denote the heavy field x, solid lines the light field (p; the shaded region 
denotes the boundary, and the dashed line the insertion of a 7-vertex. This diagram is 
easily evaluated to 



{4>(xi,yi)4>( x 2,y2)) x -effect = -2g^G N (xi, yr, x 2 , yo)<%2 - yo) 

2igj5(y 2 - yo] 



(27T) 



<Fk 



e ik x (asi -x 2 )+ik y (yi-yo)-%2 



k 2 + Ml 



. (2.23) 



Approximating the denominator in the standard way by a geometric series valid for M 2 3> A 2 , 



2igj6, 



2/2-2/0 



E 



<Fk 



2\ n 



Ml 



, 2 

ik x (xi-X2)+ik y (yi-yo)-jj 



XV "I „=o 

we extract the k y dependence in the second term as a derivative to find 17 



(2.24) 



2ig~/6. 



2/2-2/0 



M2(2vr)4 



2ig^fS. 



2/2-2/0 



E 

n=0 
00 

E 

n=0 



n 

Ml 



.a 9 a2 (*i-x2) 2 a2 (vi-vq) 2 

A 4 7r e 4 4 



(2.25) 



Now recall that the projection onto the boundary of bulk terms appears as a distribution 
with resolution A. In this sense the above term contains the delta function -^=e~^ 2< " y ~ y(l ^ !//4 . 
Up to this resolution the above expression is thus equivalent to 

2ig-fSy 2 - yo <v 



M x 

X 71=0 



E 



/ □ \ n 



(2.26) 



17 Note that these results are not inconsistent with our earlier calculation ( |2.1C| ). There we evaluate the 
answer in the approximation A S> m. Here we approximate A <C M x . The exact intermediate answer 



obtained in eq. (2.10) is non-perturbative in A/M. This is why we approximate the momentum integral 
for M x S> A in the standard way. 
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Hence we see explicitly the resultant higher derivative boundary interactions in the 
energy effective action. The above results correspond to the vertices 



low- 



s', 



eff 



d3x w & 

X 



M +0({d/M) 4 



(2.27) 



This supports the naive integrating out of x after a shift x ~^ X ~ <?(D + M 2 )~ l (f) as argued 
in section LI . The terms arising from the boundary term <f 7X<A under this shift precisely 
reproduce the higher derivative terms ( |2.27 ). 

Note the similarity between the expression ( gjgp and the image-charge term in the 



seagull-graph (2.9). We see therefore that a similar set of higher derivative corrections can 
arise from /oop-diagrams in a x<fi theory with only the bulk interaction 



qint 
°high 



tftxdy 



— 2/2 
9X <P 



(2.28) 



This is the hybrid inflation inspired model, considered before in the context of decoupling in 
FRW-spacetimes [f7[ . The seagull diagram responsible for the higher-derivative corrections 



is a direct copy of eq. (2.9) only to be evaluated in the limit M x 3> A rather than ms <C A. 




(4>(x 1 ,yi)(p(x2,y2))x-effect 



-igG N (x 1 ,y 1 ;x 1 ,y 1 )5 3 (x 1 - x 2 )S(y 1 - y 2 ) 

-fl'€;l,2 ( W 



(2tt)' 



d A k 



e a 2 



k 2 + M2 



d kx dky ~ 



Jk y (-2y+2y )-j^ 

k% + kl + M 2 



(2.29) 



Repeating the geometric series expansion in k 2 /M 2 



-g^ ; i,2^y;i,2 
M 2 (2iry 



E 

n=0 



, —h 
d 4 k' 



2\ n 



e a* + 



J d k%dky 



h z — h 2 



Jk y (-2y+2y )-fj 



(2.30) 



we see that we can extract the k y dependence in the second term as a derivative. The 
x dependence along the boundary and the full bulk term give purely local corrections as 
expected from loop graphs. Though this non-local y-dependence is counterintuitive, the 
physical reason is easily identified. It is the interaction with the image charge. We find 



bulk + 



bulk + 



1.2 



M|(2vr)4 

-g^;l,2^;l,2A 3 

M2(2vr) 4 



EEL) ( M2 



d 



2 \P 



d dk' 



x u,rv y 



h 2 \ ^ f,2 



Ml 



EE 

n=0 p=0 



n—p 



n 



P 



a: 



M 2 X 



" 1 1 dkj k * { - 2y+2m) ~^ 
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bulk + 



12 tf„;l, 2 AV/2 



M2(2vr) 4 



EE« 



n—p 




3-Y Ae -^ 2 (y-yo) 2 



(2.31) 



n=0 p=0 

where a n = 27r 3//2 (— 2) n+1 (2n + 1)!!. In the distributional sense this is therefore equal to 



bulk + 



p=0 Ivl x 



d 2p 
V S(y 



yo) 



(2.32) 



where Q p can be read off from ( 2.31 ). The bulk one-loop ^-diagrams therefore gives rise to 
the higher-derivative irrelevant corrections on the boundary 



S, 



eff 



E 



d 6 x 



g{3 p A 3 1 ( dl p 



Ml 



(2.33) 



This result shows that the boundary irrelevant operators will generically not appear in the 
combination § di4>di4> — 4>d 2 (p. This is a direct consequence of the fact that the boundary 
breaks Lorentz invariance. Examples which generate the other two irrelevant operators are 
easily found. The model just discussed will also generate § 4> A terms. A non-linear sigma 
model will naturally have § d n (j)d n (j) corrections. 

2.4.1 Minkowski space boundary conditions and irrelevant operators 

An important question therefore is how generic the occurrence of irrelevant corrections is. 
In particular do fixed points of boundary RG-flow, e.g. the Minkowski boundary conditions 
or other 'vacua', still receive irrelevant corrections? RG principles tell us that we should 
expect them. Just because we are at a fixed point of RG-flow, does not mean that irrelevant 
operators encoding a high-energy sector are forbidden. In the context of boundary RG-flow, 
the connection between boundary conditions and 'vacua', makes this statement somewhat 
surprising. In Minkowski space in particular we do not expect that integrating out a high- 
energy sector would change the vacuum state in the low-energy effective theory even at the 
irrelevant level. 18 Both the general RG principles and the intuition that in Minkowski space 
high energy physics should not change the low-energy boundary conditions are true, as we 
will now illustrate. The first point is evident from the two scalar theory at the beginning 
of this section with the interactions given in ( |2.21 ). Integrating out the % field exactly, 
clearly gives rise to the following irrelevant contributions to the low-energy effective theory 
for 6. 



qint 
low— energy 



d 3 xdy - <j)(g + j5(y - y ))(\J bc 



Ml 



n=0 J X \ 



□ 



Ml 



+ 



Ml 



1 (g + i$(y - vq))4> 
n 6cv 



Ml 



5(0)0 . (2.34) 



Here n& Cy should be interpreted as acting on a complete set of eigenfunctions with the 



boundary conditions d n x 



-kx that belong to the massive field x- To address the 



3 We thank Jim Cline for emphasizing this point. 
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formal divergence of the delta function at its origin, 6(0), recall first that in a cut-off 
theory, as we are considering, all distributions become smeared on the scale of the cut- 
off. The 5(0) in the second term is therefore proportional to M x purely on dimensional 
grounds. Our cut-off scheme eq. ( |2.10 ) indicates that S(x) = liniA^oo 7r _1 / 2 Ae~ A " ' ~ 



8(0) = Mit~ 1 / 2 . This regularization only postpones the problem, however. In appendix [D| 
we perform a computation, which indicates that the 5(0) term arising from discontinuous 
field redefinitions does not explicitly appear in bulk correlation funcions. Its sole function is 
to generalize all distributions so that they obey the correct boundary conditions d n f(y) = 
-Kf(y). 

Consistent with the principles of decoupling, we see that whatever boundary conditions 
we choose for <fi including fixed points of RG flow, the boundary action will receive irrelevant 
corrections. How can this possibly square with the idea that Minkowski space high energy 
physics should not correct the vacuum choice, i.e. the Minkowski space boundary conditions 
of <pl In this simple model it is fairly easy to see that the boundary conditions of <fi change, 
because the massive field x does not have Minkowski space boundary conditions. When 
X is integrated out, this reverberates in the low energy effective boundary action for (p. A 
naive way to see that x is n °t a t a fixed point of boundary RG-flow, is to note that the full 
boundary condition for x reads d n x = —kx — 7</>- The explicit dependence on 4> perturbs 
one away from a x-sector fixed point Kfi xe d- To consider a fixed point in the x-sector alone 
is inconsistent of course; the full x~§ dynamics needs to be taken into account. But an 
exact answer, possible because the theory is exactly solvable, shows that this naive guess 
is qualitatively correct. The exact answer is obtained by diagonalizing the theory to two 
fields $1 and $2 with action 

Smh = \j d?xdy *! (u - Ml + J^j $ a + d> 2 (u - ml + Jj^j d> 2 + 0(g*) , 

$ 2 + 0( 1 \g 1 2 ,g 2 1 ) , 



Shound = ~ <f> d 6 X $1 -^4- + 1 V 2 $1 - $2 -^T + 

bound 2 J \Ml Ml J \Ml Ml J 



Ml = Ml- ml. (2.35) 

If we tune 7 and g such that one of the two fields has Minkowski boundary conditions 
= — iu(k, M§ 2 ), we see that the difference in masses M$ 1 ~ M x and M§ 2 ~ m^ 
prevents the other from obeying Minkowski boundary conditions. 

At a very fundamental level these results are easily understood. Recall that the 
Minkowski boundary conditions are the only boundary conditions respecting Lorentz in- 
variance; this is what guarantees that the values of the boundary couplings correspond 
to a fixed point. The explicit boundary interaction § —JX'P — ~\ J S(y ~ Uo^X'fi breaks 
Lorentz invariance, however. In the diagonal system with $1, < I>2, the Lorentz invariance 
is broken because one of the two fields does not obey Minkowski boundary conditions. 

We have only shown that irrelevant operators will generically appear in a situation 
where a field in the high energy sector is not in the Minkowski vacuum. Lorentz symmetry 
should guarantee the converse: that if all massive fields obey Minkowski boundary con- 
ditions, no boundary RG-flow or boundary irrelevant operators can appear. Importantly, 
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in the setting of interest to us, FRW cosmology, Lorentz invariance is absent. It is there- 
fore not clear that cosmological boundary conditions, to which we turn now, are similarly 
protected from RG-flow and irrelevant contributions from high energy physics. Strictly 
applying the RG principles, we should not expect them to be protected. 

3. Boundary conditions in cosmological effective Lagrangians 

We have seen that: 

(1) a boundary action can encode the boundary conditions one wishes to impose on the 
fields. 

(2) This holds in full generality. The boundary need not correspond to a physical ob- 
struction or object. Completely transparent boundary conditions exist that mimick 
the situation as if there is no boundary. Introducing a boundary action to account 
for initial conditions therefore places no additional constraints on the theory. 

(3) Generically the boundary conditions will be affected by RG flow, and suffer irrelevant 
corrections that are controlled by the high energy physics. 

We now use this knowledge to describe FRW cosmologies from a Lagrangian point of view. 
The main issue in the Hamiltonian description of FRW cosmologies is that of vacuum 
selection. In the absence of a global time-like Killing vector or asymptotic flatness, there 
is no unique vacuum state. There are two preferred candidates, the Bunch-Davies and 
the set of adiabatic vacuum states, which we review below, but some uncertainty remains. 
Whichever state is the true one, points (1) and (2) above tell us that we can account for 
this state by the introduction of a specific boundary condition at an arbitrary time to- 

Our lack of knowledge of the specifics of the very early universe and the high energy 
degrees of freedom dominating at that time rather suggests to encode the initial state 
uncertainty in a 'past boundary' for any cosmological theory. With the boundary comes 
the Lagrangian translation of the vacuum choice ambiguity: what boundary conditions to 
impose? We will not give an answer to this long-standing question. We will show, however, 
that whatever (local relevant) boundary conditions one chooses, they are consistent in the 
sense that the backreaction is under control. The counterterms appropriate to the boundary 
conditions specified that are necessary to render the Minkowski stress-tensor finite, do so 
in cosmological setting as well. This confirms the intuition that the boundary conditions 
do not affect UV-physics. And this continues to hold for any choice of cosmological initial 
conditions. This may come as a surprise. The Hadamard condition — that at short dis- 
tances the two-point correlation function is the appropriate power of the geodesic distance 
o~(x\, X2) d ~ 2 — has long been thought to be a consistency requirement for cosmological 
boundary conditions. Only these correlation functions permit 'renormalization' by the 
standard Minkowski stress tensor. The lesson from section [|, however, is that other short 
distance behavior does not necessarily signal an inconsistency, but instead implies that the 
'boundary conditions' need to be renormalized as well. This returns to the front the ques- 
tion which boundary conditions describe the physics of the real world, but none that can 
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be deduced from local relevant boundary interactions are intrinsically inconsistent. This is 
the power of the effective Lagrangian point of view. 

Suppose for now that all choices for boundary conditions on the initial surface of 
an FRW universe are indeed consistent. Compared to Minkowski spacetime there is a 
new ingredient. The boundary condition needs to be covariantized. This is done by the 
introduction of a unit vector normal to the boundary. 

d n <j) = n»d^ = , \g^n u \ = 1. (3.1) 

In the conformal frame, 

ds 2 FRW = a 2 (i]){-di] 2 + dx 2 d _ x ), (3.2) 
the unit normal vector to the boundary scales as a . Hence the boundary condition reads 

~9ri4>\ri=rio = —K4>\r]=r] • (3-3) 

The explicit dependence on the scale factor a simply reflects that momenta redshift under 
cosmic expansion. 19 To construct the two-point correlation function for a massive scalar 
that satisfies this boundary condition, we need the equation of motion in an FRW back- 
ground. For simplicity we will assume that this background is pure de Sitter; the results 
below generalize straightforwardly to power-law inflation and are therefore truly generic. 
The equation of motion is 



1 



dfiV^gg^du^x, 77) - m (f>(x, rj) = 



V~9 

( 1 a i p \ 

=> -^dl + (d - 2W9„ + -= + m 2 ) <f>(k, rj) = . (3.4) 

In the second step we Fourier transformed the spatial directions. Substituting the constant 
de Sitter Hubble radius a~ 2 a' = H, the explicit scale factor a = —1/Hrj and making the 
conventional redefinition rj = —y/k, we have a Bessel equation for <j) = y^^ 1 ^ 2 ^'- 

fd 2 + ydy + y 2 + ^- M y) = 0. (3.5) 

The most general solution to the field equation is therefore 

W s, + S (-k^-"'\fZ(")^HA-h). „ = ,/00£, (3.6) 



ik\JJ ' V 4 IP 



with H v {y) the Hankel function satisfying eq.( |3.5|) . The normalization and convention is 
such that in the limit k — > 00 we recover the Minkowski space solutions. The boundary 
conditions (|Q|) determine b, as in eq. (|2,13|). 



19 Realizing that cosmological scaling induces RG-flow we manifestly see the previous claim that Dirichlet 
conditions are trivial IR-fixed points. 
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By construction the Green's function is given by (see appendix [A| for details) 20 

G KfiK {k\,r]x]k2,r]2) = (2vr) 3 5 3 (A;i + k\)N Kf , K \^b Kf (ki,m) < Pb K (^ T t2)0(rii ~ V2) 

+ip K (k 1 ,r] 1 )Tp b Xk2,mW(V2-Vi)) , (3-7) 



where Kf characterizes the future boundary conditions at y = 00. The normalization M K . K 
is chosen such that (□ — m 2 )G = i5 d /^—g. This requires that 

M Kf , K n K (k, r]) d v lp hR! (k, 77) = -ia 2 ~ d ( V ) = -i(-Hr ] ) d - 2 . (3.8) 

We find that 

1 

(1 -J Kf b K ) 

From here on we will again restrict our attention to d = 4 spacetime dimensions. 

3.1 Harmonic oscillator and shortest length boundary conditions 

A special set of boundary conditions are the covariantization of the completely transparent 
"Minkowski" boundary conditions of eq. ( 2.15| ). We will call these "harmonic oscillator" 



M Kf:K = =^-^ . (3.9) 



boundary conditions. Recall that these correspond to the boundary action § (j)J df — m 2 (f>. 
Covariance requires that the scale factor should enter here as well. We thus find that the 
cosmological harmonic oscillator boundary condition is characterized by 



k 2 

kho = -f +m 2 . (3.10) 



a, 







For the specific momentum dependent choice of boundary location TjQ L (k) = — A/H\k\ 
or equivalently ao = \k\/A, these boundary conditions correspond to a constant value 
for the physical parameter b. They are therefore the boundary conditions proposed in 
J2|, Underlying this inspired choice is the thought that in a cosmological theory there 
is an 'earliest time', where a physical momentum p = k/a(r)) reaches the cut-off scale (the 
shortest length). Whether there is truly an earliest time in cosmological theories is an 
interesting question in its own right. It would be the natural location for the boundary 



A 'covariant' Green's function is given by 

trunc(n j ) 



4>b K ,n(ril)4>b K ,n(ri2) 

U) H 2 n 2 - m 2 + H 2 (d - l) 2 /4 



where Kf characterizes the future boundary condition at rj — 00 and fJ,(n) is an easily determined measure. 
From this expression it is clear that the delta function therefore also obeys the boundary condition. Indeed 
the delta function is best viewed as a completeness relation for eigenfunctions of the Laplacian = —k 2 (p 
obeying a^ 1 d n (p k \ Vo = -K(fk\ vo , i-e. 
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action, but as a consequence of the symmetry between boundary location r/o and coupling 



k exposed in section 2.2, it is not directly relevant to us. Indeed it is easy to see that a 
momentum-independent coupling kho at i]Q L (k) = —A/H\k\ is equivalent to a boundary 
action on a standard time-slice rj' with momentum-dependent coupling ksl 

KSL ~ cj> + (ri'o) + bsL<i>-(v'a) ' SL " + ^-(^) ' ^ 

In the limit A — > oo we recover the harmonic oscillator vacuum at rj = — oo. The 
coupling k' encodes these harmonic oscillator boundary conditions at rjo = — oo in terms of 
conditions at rj' plus corrections that vanish as A —* oo. As we have seen in the previous 
section and will discuss in detail in the next, these corrections therefore correspond to the 
introduction of specific irrelevant boundary operators. 

3.2 The Bunch-Davies and adiabatic boundary conditions 

In universes without a global timelike Killing vector, there is no clear concept of the 
vacuum as a lowest energy state. Particle number is also not conserved and one cannot 
unambiguously define an 'empty' state either. Instead one must specify a particular in-state 
characterizing the initial conditions. Two solutions to this vacuum choice ambiguity have 
become preferred. One is the Bunch-Davies vacuum, which is indirectly constructed by 
requiring that for high momenta k/a^> H the Green's function reduces to the Minkowski 
one. The second corresponds to the set of (n-th order) adiabatic vacua, which is constructed 
by the requirement that the number operator on the vacuum changes as slowly as possible 
g H). For de Sitter space the infinite order vacuum and the Bunch-Davies one are the 
same; we shall therefore only discuss the latter. 

The boundary conditions corresponding to the Bunch-Davies vacuum are readily found. 
In the basis ( |3.6| ) we have chosen, the Bunch-Davies-state corresponds to choosing 6 = 0, 
and hence 

K B D = — • (3.12) 

fdS,+,0 

Note that the Bunch-Davies boundary conditions are the analogues of the Minkowski 
boundary conditions in a mathematical sense only. The flat space Minkowski boundary 
conditions in eq. (j2.16|) are easily recognized as K M f~ k space = -d n f Mink, +,o/<p 'Mink, +,o with 
<PMink,± — e^ lut . Using the Bessel function recursion relation 

d y H v (y) = -H u (y) - H u+1 , (3.13) 

and the chain rule d v = —kd y (recall that d n = a~ 1 d v ) a straightforward calculation yields 

, <r+i(-kvo) , (d-l) + 2v\ 
kbd = — , — - — ; h 




H v (-krjo) 2krj 
Iv+i{-kr]Q 
H v (-kr) ) 



H u+1 (-k Vo ) \ + H (d-l) + 2v ^ (314) 
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Knowing the asymptotes of the Hankel functions 

1 (lY r(i/) (2\ u , . 

z : H u {z) - _ , 3.15 

sm(z/7r)r(l — i/J \Z/ IT \z J 

z^oo: H u (z) ~ J —e^-h^-z*) , (3.16) 

V 7TZ 

we see that for 770 — > — 00 the Bunch-Davies boundary condition reduces to harmonic 
oscillator boundary conditions 

|fc| / i*\ , „(d-l) + 2i/ 

~ -i^ (3.17) 
a 

of a massless field. (One cannot say that the boundary conditions tend to Dirichlet, the 
diverging ao is compensated by the normal vector, see eq. (|3.3|) .) The mass correction is 
subleading in this limit. We should keep in mind though that this is a formal expression. 
At rjo = —00 the induced boundary volume vanishes, and boundary conditions cannot 
easily be accounted for in terms of a boundary action. 

3.3 Transparent, thermal, adiabatic boundary conditions; fixed points of bound- 
ary RG flow? 

The most natural choice for the boundary conditions are arguably the ones which are 
transparent. If there is no real interface at the boundary location yQ, no physical effects 
of its location should be noticeable. To define transparency we need a notion of incoming 
and outgoing waves. A clean definition of such waves only exists in asymptotically flat 
spaces. Suppose one establishes these and let us call the incoming wave (from the past) 
and the outgoing 92+. The transparent boundary conditions are then those with b K = 0. 
Of course de Sitter space is not asymptotically flat, but based on the asymptotic behavior 
of the Bessel functions, one can argue that the basis functions <fds,- an d Vds , ,+ defined in 
Q3.6D correspond to in- and out-going waves respectively. In that sense the Bunch-Davies 
boundary conditions are the transparent ones. 

A definition which is more intrinsic to de Sitter is that the Bunch-Davies boundary 
conditions are the thermal boundary conditions. This emphasizes the existence of a cosmo- 
logical horizon, and is probably tied to the notion of transparency. From the Lagrangian 
point of view the true vacuum should be a (UV) fixed point of boundary RG-flow. In the 
presence of a global timelike Killing vector with a conserved quantum number dt<j) = iE<p 
such a fixed point is easily constructed following the Minkowski space example in section 



2.3. In cosmological spacetimes it is not clear what the fixed points of boundary RG-flow 
are or whether there are any. The absence of a unique vacuum suggests that there may 
be none. If we recall that cosmological expansion induces RG-flow, the definition of the 
adiabatic vacuum, i.e. that the number operator on the vacuum change as slowly as pos- 
sible, becomes very interesting. It would be worthwhile to investigate these connections 
between the transparent (i.e. Bunch-Davies), the thermal, and the adiabatic vacuum in 
FRW backgrounds and fixed points of boundary RG-flow further. 



-21- 



3.4 Backreaction and renormalizability for arbitrary boundary conditions 

We shall now make a crucial point. Any cosmological boundary condition n, provided 
it is a dimension-one analytic function of the spatial momenta, is consistent in the sense 
that backreaction is under control. The divergences appearing in the stress tensor must, 
of course, be regulated by the flat space counterterms of the same theory. This includes 
the boundary counterterms for <f kc/) 2 and <f fj,(j)d n (fr. Our review in section ||| has made 
this clear. In a rather coarse fashion we can also see this directly from the FRW Green's 
function in the limit of high (spatial) momentum — in as far as this limit exists in a cut-off 
theory. Using the asymptotic values of the Hankel functions, the basis functions 4>±,ds{k, T}) 
tend to massless Minkowski ones (the mass is negligible in the high momentum limit) 

t a ft \ 1 e±ikV <f>±,Mink(k,r]) 

k -> oo : <p±,ds{k,ri) ~ — = = — ■ . (3.18) 

v 2k ^ ^ 

The coefficient b encoding the effective boundary conditions for high-momentum modes 
therefore does not vanish, but reads 

, K(j}+ t Mink,0 + a o 1 d v 4'+,Mink,0 ~ H(f) +j Mink,0 



Mink,0 + a o ^0-,Mmfe,O — H(j)_Minkfi 

e 2i\k\r l0 _ (3 19 ) 



a K + i\k\ + a H o2ilklvo 



ciqk — i\k\ + aoH 

The last terms in the numerator and the denominator are negligible in this limit \k\ 3> 
aH. They are remnants of the fact that the background breaks Lorentz invariance. The 
coefficient b thus does not vanish in the high momentum limit. Because a non-zero b means 
that there will be divergences in the theory aside from the 'Minkowski '-space divergences, 
it appears that any choice of boundary conditions with b ^ is in trouble. In section || 
we reviewed, however, that this is not so. The additional divergences are localized on the 
boundary surface where the boundary conditions are imposed, and can be reabsorbed in 
a redefinition of the boundary couplings. Any choice for b (descending from a boundary 
coupling k that is dimension one and analytic in the spatial momenta) is consistent. 

One is tempted to conclude that for any boundary condition imposed at rjo = — oo, the 
high spatial momentum limit of b vanishes. This is true in the sense that if we keep k fixed 
our flat space intuition, that boundary effects vanish when the boundary is moved off to 
infinity, continues to hold. However, this goes against the principles behind the framework 
we advocate here. In the sense of the symmetry between boundary location and boundary 
coupling k, as explained in section [T^, it is only the specific combination b K which matters. 
At what location r/Q one imposes the boundary conditions n is immaterial to the physics. 

The conclusion is that the answer to the question "what boundary conditions should 
we impose on quantum fields in FRW backgrounds" requires physics input rather than 
internal consistency. The Bunch-Davies vacuum certainly seems the closest analogue of 
Minkowski boundary conditions, even though it is not the naive covariantization of them. 
The similarity suggests that the Bunch-Davies boundary conditions may correspond to a 
fixed point of boundary RG-flow. At the same time Lorentz symmetry is still broken. If 
they are renormalized, it would suggest that they are not special in any sense. 



- 22 - 



4. Transplanckian effects in Inflation 



Inflationary cosmologies are the leading candidates to solve the horizon and flatness prob- 
lems of the Standard Model of Cosmology. Consistency with the observed spectrum of 
temperature fluctuations in the Cosmic Microwave Background (CMB) provides an esti- 
mate of the Hubble parameter H during inflation. Depending on the model, H can be as 
high as 10 14 GeV. With the string scale M s t r ing = 10 16 GeV as the scale of new physics, 
this means that the suppression factor H/M of irrelevant operators could optimistically 
be at the one-percent level. This opens a window of opportunity to experimentally witness 
effects of Planck scale physics [|J. Besides its theoretical appeal, inflation is also the leading 
candidate for early universe cosmology on experimental grounds. The most precise cosmo- 
logical measurements to date, the temperature fluctuations in the CMB, advocate inflation. 
The CMB measurements are therefore also the most promising arena where remnants of 
transplanckian physics could show up. In inflationary cosmologies the CMB temperature 
fluctuations originate in quantum fluctuations during the inflationary era. The issue of 
vacuum selection in cosmological settings thus has immediate consequences for CMB pre- 
dictions. At the classical level the Bunch-Davies choice is, for reasons reviewed in the 
previous section, the preferred one; it is the closest analogue to the Minkowski boundary 
conditions. Previous investigations into effects of Planck scale physics suggest that the 
CMB fluctuation spectrum is affected at leading order in H/Mpi anc k and that this effect is 
precisely due to the choice of vacuum [||, || . Due to our ignorance of the details of Planck 
scale physics (i.e. our lack of understanding of string theory in time-dependent settings), 
decoupling in effective field theory is arguably the framework in which transplanckian cor- 
rections must ultimately be understood Q. By the addition of an arbitrary boundary 
action encoding the boundary conditions, we have put the issue of vacuum selection on a 
consistent footing with the ideas of effective field theory. In this comprehensive formula- 
tion, we can deduce systematically what the effect of Planck scale physics is on boundary 
conditions (vacuum selection) and whether its effect on CMB predictions is indeed leading 
compared to bulk corrections. 21 

The Planck scale physics is encoded in irrelevant operators. The leading bulk irrelevant 
operator J cj)D 2 (j) consistent with the symmetries is dimension six. In section 2A we 
constructed and derived the four leading irrelevant boundary operators in flat space 

' d 3 x ^ , — / d 3 x d^dicj) , — / (fix d n (^d n (t> , — / d 3 x 4>d n d n $ . (4.1) 



Ml Ml Ml M , 

J y=yo 1V± J y=yo 2,1 J y=yo 1V1 J y=yo 

compatible with unbroken I SO (3) symmetry. In a cosmological setting this is the require- 
ment of homogeneity and isotropy. These operators are all dimension four and as the 
explicit scaling shows, they are expected to be dominant over the leading bulk irrelevant 
operator. In curved space these operators are covariantized. For a scalar field (ft covarianti- 
zation has only a significant effect on the last operator in (14.11). A new coupling is needed 



21 The object of our study is an external scalar field in a fixed FRW background. Strictly speaking only 
the gravitational tensor fluctuations are effectively described by such a model. However, our arguments 
should apply to the scalar-metric fluctuations as well, since these only differ by an amplification factor of 
the inverse slow-roll parameter. 
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which provides the connection for the covariant normal derivative 

^ j sfhn»n v {^d^ - W^dpigW) = j> VhnWD^ . (4.2) 

Here hij = g uu diX fl djX u is the induced metric on the boundary, and its unit normal 
vector. In FRW cosmology with the metric in the conformal gauge, 

ds 2 FRW = a 2 ( V )(-dr] 2 + dxl) , (4.3) 

and an initial timeslice r/ = r/Q as boundary, the induced metric, connection coefficients, 
and normal vector are 



h 



iv 



a (S, 



y ) 1 

°0 1 



a H 5. 



T l ■ 



(4.4) 



Here clq = a(%) and Hq = H(t]q) is the Hubble radius H = a 2 d r) a at r/ = i]q. Substituting 
these values we obtain the FRW version of the irrelevant operator 



j alcf) (d n - H) d n(l 



(4.5) 



We shall compute the effect of the leading irrelevant operators on the two-point corre- 
lator of (p. In inflationary cosmologies, the latter determines the power spectrum of CMB 
density perturbations. We will assume we can treat the four-point bulk Xcj) and (irrelevant) 
boundary interaction <f cj) 4 perturbatively and will ignore them to first order. Combining 
the remaining irrelevant boundary operators in a correction to the FRW boundary action, 
one obtains 



qirr.op. 
bound 



V=V0 



2M 



d' 



(4.6) 



The precise value of a coupling constants (3{ is determined by two parts. (1) It is determined 
by the details of the transplanckian physics; e.g. if transplanckian physics is a free sector, 
decoupling is exact and (3 = (for dynamical gravity the sectors are never decoupled of 
course), but (2) the couplings f3{ are also covariant under the symmetry between boundary 
location and coupling. If we would have computed the irrelevant corrections to a boundary 
condition at a different location i/q, we would have found different values which upheld 
that all physical quantities only depend on the choice of boundary location through a 
specific combination h K ^. 



Two of the operators in eq. (^6) contain normal derivatives. As discussed in section 
|2[ such operators can be removed by a discontinuous field redefinition and a change of the 
remaining boundary couplings. We do so in appendix |C]. To lowest order in 0i/M, eq. 
(|4.6|) is equivalent to a boundary interaction (if the boundary coupling /U=0) 



girr,leading 



a 3 d 3 x - 



+ 



M 



M 



(4.7) 
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where m 2 is the mass of the scalar field. Fourier transforming along the boundary, the 
leading irrelevant correction thus amounts to a change in the boundary condition k by 22 

K ^f = ^+ alM +^ r -^ r -^^ r - ( 4 - 8 ) 

We clearly see that the leading correction to the low-energy effective action occurs at order 
\k\/aoM and H/M. For CMB physics the momentum scale of interest is \k\/a presen t ~ H, 
and both are of the same order. The conclusion that the \k\ dependent operators are 
suppressed by a factor ao/a presen t is incorrect, when we recall that the location of the 
boundary is arbitrary. 

For a given FRW universe the Green's function, including the H/M correction to 
the boundary condition, can now simply be read off from eqs. ( |3,6| )-( |3~77| ), We can thus 
straightforwardly compute the leading transplanckian effect on the power spectrum of 
inflationary perturbations. The latter is related to the equal time Green's function with 
Kf = k (see appendix [E|) 

k 3 - 
P(k) K = lim ——G Kf= ^ K (k, n; —k, rj) 

= l im j M^(M)! 2 (49) 
J™ 2tt 2 (1 - |6„|2) ' l4 ' 9j 

where (/c, 77) is a solution to the (free) equation of motion, normalized according to 
the inner product (3.8), and with boundary condition d n f \ = — Ktp\. Note that the basis 



functions (fb K only depend on the location of the boundary through the physical combination 
b K . This 'independence' of the location of the boundary guarantees that the power- spectrum 
- a physical quantity — is so as well. For an infinitesimal change in the boundary 
condition k, we can treat the vertex § — ^Skc/) 2 perturbatively, and the change in the power 
spectrum simply amounts to computing the following Feynman diagram. 




(4.10) 

This immediately illustrates that if 5k is of order H/M, the change in the power spectrum 
will be of order H/M. For completeness, we compute the power spectrum by de Sitter 
Feynman diagrams in appendix ^. With the effective change in k corresponding to the 
contributions of the irrelevant operators f3{ known, we can also simply expand the exact 

22 Because the coupling k is subject to renormalization, its value is fixed by a renormalization condition 
and an experimental measurement. An important question therefore is, whether the effects of irrelevant 
operators are experimentally measurable. The standard story, that (1) measured couplings always include 
all relevant and irrelevant corrections, and that (2) the contribution of each coupling f3i is an independent 
contribution to the precise running of coupling n e ff(f3i) under RG-flow, should apply. A very precise 
measurement of the scaling behaviour of k should reveal the contributions of high energy physics encoded 
in the irrelevant operators. 
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solution for the power spectrum for any n. Choosing the Hankel functions as basis as in 
eq. ( |3,6|) , the solutions iff, K are given by 



For an infinitesimal shift 5k the power spectrum is thus 



(4.11) 



k 

P(k) K+SK = P(k) K + lim — 
v^o 2ir z 



5b 



FW^ + c ' c ' 



+ 0(5b 2 



(4.12) 



Substituting the de Sitter values computed in the previous section, and using that asymp- 
totically (see ( 3. 15| ) ) 



r - v 



(4.13) 



we find that 



P(k) K+SK = P K 1 + 



(i-H 2 ) 2 



*&4 + a t 



(4.14) 



Recall from eq. ( |2.14 ) that 



5b 



+ 



(4.15) 



We see explicitly that the change in the power spectrum is also linear in H/M. 

For the preferred Bunch-Davies vacuum choice, where 6 = 0, the corrections thus 
become 



P B D+5n{k) =Pbd\1 + 



6k- 



-<t>-,od n <t>+,o + (f>+,od n <P-,o 



+ c.c. 



(4.16) 



It appears we have introduced a dependence on the boundary location, but we should 
not forget that 5k intrinsically depends on yo as well. The combination above is guaran- 
teed to be independent of the boundary location. We recognize in the denominator the 
normalization condition ( |3,8| ) (with d n = a~ l d v ). The expression therefore simplifies to 



Pbd+Sk = Pbd\1 + 



5k- 



c.c. 



-ia r 



+ 0(5k 2 ) 



(4.17) 



Restricting our attention to de Sitter space, we insert the explicit expressions for the basis 
functions 4> + from eq. ( |3.6|) , and obtain, using that oq = k/Hyo, 



rydS T)dS I i 

BD+Sk — ^BD 1 



(w) 



6KH v (y ) 



+ c.c 



(4.18) 
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Substituting the irrelevant operator induced 5k from eq. ( [4.8[) , we compute the following 
corrections to the power spectrum 



pdS _ pdS 

^BD+Sk r BD 



7T 

w 



Hl(y ) 



-/5c) , k 2 bd P ± f3 c m 2 3[3 C H 

H — — K BD 



M 



M 



M 



+ c.c. 



(4.19) 



with (eq. flO^ )) 



K BD 



d-l + 2v 



H 



k H 



v+l 



(w>) 



(4.20) 



This is our final result. Let us stress again, that the apparent dependence on the bound- 
ary location is only that. The boundary couplings by construction compensate the y$ 
dependence and the whole expression is independent of yo- 



5. Conclusion and Outlook 

The recent successes in CMB measurements exemplified by || , have made the computation 
of inflationary density perturbations a focal point of research. The computation of these 
density perturbations suffers from a fundamental deficiency, however, that is at the same 
time a wondrous opportunity. The enormous cosmological redshifts push the energy levels 
beyond the bound of validity of general relativity, the framework in which these compu- 
tations are done. From a field theoretic point of view general relativity can be viewed 
as the low energy effective action of a more fundamental consistent theory of quantum 
gravity. This effective action has higher order corrections which when re-included increase 
its range of validity. These higher order corrections encode the physics that is specific to 
quantum gravity. Hence understanding the way these higher order corrections affect the 
computation of inflationary density perturbations is both needed to restore consistency to 
the computation, and provides an opportunity to witness glimpses of Planck scale physics 
in a measurable quantity. 

However, an action by itself is not sufficient to extract the physics of quantum fields. 
One must in addition specify a set of boundary conditions. Which boundary conditions to 
impose is always a physical question. In the Hamiltonian language boundary conditions 
correspond to a choice of vacuum state. In cosmological settings, due to the lack of sym- 
metries the correct choice of vacuum, i.e. boundary conditions, is ambiguous. A number 
of proposals, though, exist for the correct state. What we have discussed here, is that this 
vacuum choice ambiguity can be framed in terms of the arbitrariness of a boundary action. 
This puts the full physics in the form of a naturally coherent effective action. Deriving 
the power spectrum of inflationary density perturbations within this framework, the lowest 
order corrections are irrelevant boundary operators of order H jMpi anc ^. Because we are 
able to use the language of effective field theory, not only is the parametric dependence 
of the inflationary perturbation spectrum on high-energy physics known, the coefficients 
are also in principle computable from the high-energy sector that has been integrated out. 
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RG-principles tell us that generically this coefficient will be non-zero, except for very spe- 
cial choices of initial conditions and high energy completions of the low energy theory. In 
cosmological spacetimes in particular the Lorentz symmetry which forbids the appearance 
of such corrections in flat Minkowski space is absent. This makes the prediction that we can 
potentially observe Planck scale physics in the cosmic sky quite strong, or equivalently the 
absence of these effects would constrain the possible high energy completions, i.e. string 
theory. 

Several earlier investigations have shown that the effects related to a choice of initial 
conditions are not the only way in which high-energy physics can show up in cosmological 
measurements. Effects due to a non-vanishing classical expectation value of high- @ or 
low-energy || fields, or a modified dispersion relation (see, e.g. JlJ) can be of the same 
order. The former two should fit into our framework by the explicit introduction of sources. 
The latter presumes an all-order effective action, which is finite and therefore has a specific 
kinetic term J-(U\/K). The subleading effects in A obviously change the two-point corre- 
lation function and hence the power spectrum. In RG-terms a specific choice of regulator 
function ^"(D/A) corresponds to a specific choice of UV-completion of the theory. The rel- 
evant behaviour is universal and independent of the choice of ■?-"(□ /A), but the irrelevant 
corrections are not, of course. 

The introduction of a boundary action to account for the initial conditions, and its 
behaviour under RG-flow including irrelevant corrections begs for a comparison with the 
idea of holography. The latter suggests that (gravitational) theories in d-dimensional de 
Sitter space have a dual formulation as a (Euclidean boundary) conformal field theory of 
dimension d — 1 [27, 28]. The cosmological implications of this conjectured correspondence 
underline the universality and robustness of predictions for inflationary density perturba- 
tions precisely because they are related to RG characteristics in the dual d — 1 dimensional 
theory pOj , [29| , [30[| . These qualitative similarities are striking, but there are crucial differ- 
ences with the approach put forth here. Holography interchanges the IR and UV properties 
of the dual theories. The UV physics of a three-dimensional Euclidean field theory corre- 
sponds to the IR of the four-dimensional de Sitter gravity and vice versa. The holographic 
screen where the dual field theory lives corresponds to a boundary action in the de Sitter 
future. Its precise position defines the UV cut-off in the Euclidean field theory that should 
completely describe the infinite interior (i.e. the past) of the de Sitter bulk gravity theory. 
Time evolution in the bulk is then interpreted as RG-flow in the boundary field theory, and 
so the IR physics in the field theory corresponds to the infinite past in the bulk. Instead 
the boundary actions considered in this paper are introduced only to encode the initial 
conditions in the past of the four dimensional de Sitter gravity theory. They are not dual 
descriptions of the bulk de Sitter theory, but are merely introduced as effective tools to 
describe the initial conditions in the bulk. Nevertheless, it would be very interesting to 
study how the results described in this paper should be interpreted from the point of view 
of a putative dual three-dimensional Euclidean field theory. 

That early times in cosmological theories are dominated by UV physics leads to a final 
open question. Do cut-off theories in a cosmological setting cease to be valid beyond an 
earliest time? Naively this is so, and that time would be a natural location for our boundary 
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action. The freedom, however, to impose initial conditions where-ever one wishes, means 
that we do not need to answer this question to address the issue of boundary conditions 



in FRW universes. This fact is made manifest in the symmetry (2.14) between boundary 
location yo an d boundary coupling k. Physics depends only on the invariant combination 
k{uo)- With the effective field theory description in mind, and the idea that 'vacua' are 
boundary RG fixed points, a truly interesting question is whether such boundary conditions 
exist, and if so, how they are related to the known cosmological vacuum choices. 

5.1 A comparison with previous results and the discussion on a-states 

Much discussion has taken place in the recent literature on the consistency of so-called 
a-states in de Sitter space [[|, |3l|. Initial investigations into the sensitivity of inflationary 
perturbations to high energy physics found that in pure de Sitter the leading H/M correc- 
tions to the power spectrum can be interpreted as choosing the harmonic oscillator vacuum 
(section |3.l|) at the naive earliest time rjo(k) = — A/H\k\ where the theory makes sense, 
rather than the Bunch-Davies choice [§, ||. Imposing such boundary conditions in pure de 
Sitter can equivalently be interpreted as selecting a non-trivial de Sitter invariant vacuum 
state called an a-state ||. Strictly speaking, the Shortest Length (SL) boundary conditions 
are only imposed on momentum modes below the cut-off scale A of the theory, and they 
are not true de Sitter a-states. Subject to this distinction, the purported inconsistency of 



a-states, particularly with respect to the decoupling of Planck scale physics [31], therefore 
would have major consequences. If a-states and other boundary conditions are all incon- 
sistent, all high-energy physics would have to be encoded in bulk irrelevant operators. This 
would put transplanckian effects in the CMB perturbation spectrum beyond observational 
reach. 

Let us put first, that our results form solid evidence for the presence of H/M effects 
affecting inflationary predictions for the CMB perturbation spectrum. As the explicit 
expression fl4.19|) we derive for the power spectrum shows, our results, though qualitatively 
similar, are quantitatively far more general from having 'chosen' an (cut-off) a-state. The 
coherent effective Lagrangian approach followed here gives a precise answer which differs in 
general from the (earliest-time) a-state approach, but upholds the qualitative validity. One 
can certainly ask to what choice of 'vacuum state' our results correspond; given the physical 
parameter b K this is straigtforward to work out. The answer may be interesting from the 
point of view of Hamiltonian dynamics, but as we have shown here, in the Lagrangian 
language of boundary conditions, any initial state which can be described by a local relevant 
boundary coupling k is consistent. There is no need to know whether a-states are consistent 
to study transplanckian corrections to inflationary perturbations. 

At the same time, vacuum choices, a-states included, do correspond to boundary 
conditions. 23 And boundary conditions should not spoil decoupling, although there will 
be new effects, as we reviewed in section ^. Taking this lesson to heart, it is hard to see 
how (earliest-time) a-states could be inconsistent. A recent article [32 1 arguing for the 



23 We are grateful to Brian Greene both for emphasizing the importance in explicitly discussing the 
consistency of a- vacua and his help in resolving the issue. 
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consistency of a-vacua does not exactly follow the approach outlined here, but is very 
much in the spirit of introducing boundary counterterms. An answer, however, is provided 



by pursuing the discussion in section 3A further. The (cut-off) a-vacua correspond to 
choosing earliest-time boundary conditions in an effective theory below scale M with the 
physical parameter b$L a constant number. The precise relation is that bsL = e a - One 
then readily derives that an a-vacuum corresponds to a boundary coupling (see eq. ( |3.11| )) 

d n <f>+(rf ) + bsiA^JjQ 
KSL ~ Mrfo)+bs^-H) ' [ j 

Recall that bsL is constant. To analyze the high spatial momentum behavior, we may 
therefore approximate the modefunctions </>±(f?o) by their Minkowski counterparts. In this 
limit the boundary coupling ksl encoding a-states becomes 

,\k\ e ^' - bsLe -i\k\v' 
\k\ -> oo, K SL ^ -i • (5.2) 

The boundary coupling ksl therefore has an infinite set of poles 

\k\ = zl((2n+l)7r + i\n(bsL)) , n£Z, (5.3) 

in the momentum plane. Clearly this boundary coupling corresponds to a non-local action. 
Cut-off a-states, i.e. shortest length boundary conditions, therefore fall outside the class 
of local relevant boundary conditions we study here. But are they inconsistent? Recall 
that the original studies pi || argue that a-vacua should encode (first order) effects of 
high-energy physics in the spectrum of inflationary density perturbations. This point of 
view therefore states that by construction the boundary coupling ksl includes the effects of 
irrelevant boundary operators. We are therefore instructed to treat the non-local nature of 
the boundary coupling ksl in the low-energy effective action in the usual way. One expands 
around the origin \k\ = in the momentum plane generating a series of higher derivative 
irrelevant boundary operators with specific leading coefficients This expansion is 

valid as long as we limit the range of our effective action to the location of the first pole 
I I = 21^7 Vl 71 + ^ n ^sl\ 2 : i-e. physical momenta are constrained to the range |po| = 

|J^| < f \\nb S L\- (Eq. ( |3T9| ) gives us b SL ^ H/2Me- 2iM ^ H - in / 2 , and we recover the 
cut-off |p| < M.) The fact that the complicated pole structure of boundary couplings of 
alpha-vacua is highly specific (they ensure that (non-cut-off) a-vacua are invariant under 
de Sitter isometries) is not to the point in this perspective. It is then also clear why a-vacua 
are not renormalizable, in particular in the sense that the bare backreaction, the divergence 
in the stress tensor, is to leading order not identical to that in Minkowski space. Irrelevant 
operators correspond to non-renormalizable terms in the action. Because the pole structure 
of the boundary coupling k reveals that a-states are correctly to be interpreted as encoding 



24 It is not completely clear that this interpretation withstands close scrutiny. Most non-local terms in the 
effective action have real poles. Here we are confronted with imaginary poles. Perhaps a-vacua correspond 
to a high-energy completion with numerous unstable particles. 
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specific contributions from irrelevant operators, any correlation function computed with 
respect to the a-vacuum, includes the contribution from these irrelevant operators. It is 
therefore expected to be non-renormalizable. Obviously this does not mean that the a- 
vacua are inconsistent. As always in effective actions one must 'neglect' any contributions 
of irrelevant operators for the purposes of renormalization. They only make sense in a 
theory with a manifest cut-off |24|]. Removing the cut-off, removes the irrelevant operators. 
Indeed the a-states proposed in |2|, |(| with bsL — H/2M are naturally in accordance with 
this precept. 

In this sense, the (cut-off) a-vacua are therefore manifestly consistent in the framework 
put forth here. They simply correspond to a specific choice of leading and higher irrelevant 
boundary operators. Whatever they are is not very interesting from the perspective of 
effective field theory. 25 A specific choice for the irrelevant operators means having chosen 
a specific form for the high-energy transplanckian completion of the theory. But what this 
physics is, is precisely the knowledge we are after. 
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A. Green's functions and boundary conditions 



By definition the (real scalar) Green's function is the inverse of the (real scalar) kinetic 
operator 



2(i 



(n x -u 2 )G(x,x') = i 



5 d {x - x') 



{U x , - lo 2 )G(x, x') , G{x, x') = G(x',x) . (A.l) 



For simplicity we reduce the spacetime to one timelike direction. 

.5{t-f) 



(f^+u; 2 )G(t,t') 



(A.2) 



25 They are /3 C = -f e ~ 2iM/H , f3 {] - /3± = f e ' 2iM/H . Expanding around small \bsz.\ = H/2M and small 
fc| <C Hao, we see that 



ksl = —I — 

do 



-2iM/H 



gj*i . 

a M 



-2iM/H 



H -2iM/H r.-^BD —2iM / H 



Comparing with (i.8) we find the coefficients Pi. 

Recall that we are using the + + +— convention; i.e. in Minkowski space 
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The Green's function is thus a solution to a inhomogeneous second order differential equa- 
tion. The solution to (A. 2) is therefore not unique; we can always add a combination of 
the two linearly independent homogeneous solutions (denoted by the superscript (f))), 

+ J)$\t) = , +^)$\t) = , ^\t) * $\t) , (A.3) 

to a Green's function and obtain another Green's function. This ambiguity is resolved by 
imposing a set of boundary conditions on the Green's function. Consistency then requires 
that the delta function appearing in eq. (|A.2[) obey these boundary conditions as well. Let 
the Green's function obey the boundary condition 

d t G(t,t% =t0 = -K G{t ,t') . (A.4) 

Acting with + kq on G(t, t') in eq. ( |A.2 ), which clearly commutes with Jp- + u> 2 , we are 
forced to conclude that 

(dt> + Ko)(^+u 2 )G(t,t%, =to = -;(d t ,+K )^=^L =to = . (A.5) 
at V-SW 



The delta function on the RHS of eq. ( |A.2| ) is therefore not the naive Dirac delta-function, 
but contains extra correction terms which contribute only on the boundary. 
Limiting our attention to boundary conditions of the type eq. (A.4), i.e. 

(d t + Kf)G Kf>K0 {t,t')\ t=tf = 0, 

(d t + K )G KfjKo (t,t% =to = 0, (A.6) 

at the two boundaries at tf and to, we solve the ambiguity in the Green's function in 
this appendix. We do so because this subset of possible boundary conditions is a very 
interesting one. It contains both the canonical Neumann and Dirichlet cases, and for Z2 
symmetric scalars in effective field theory the boundary condition derived from all relevant 
boundary interactions are of this form. 

There are two standard ways to solve for the Green's function. The Hamiltonian 
way picks the timelike direction t as a the preferred one. In multiple dimensions the 
frequency to is given by the eigenvalue of the remaining spatial component of the Laplacian 
— lo 2 ^ u (x, t) = □^<I> a; (x, t). For the timelike direction, however, the Hamiltonian Green's 
function uses as building blocks the two independent homogeneous solutions to the kinetic 
operator of eq. ( |A.3| ) multiplying stepfunctions 6{t — t') and 6{t' — t). With some insight 
the boundary conditions are readily imposed. Realize that for the past boundary condition 
only the part with 6(t' — to) will contribute; while for the future boundary condition only 
the part with 9(tf — t') will contribute. Write the Green's function as 

G*™M = Mi (^(t)^{t')6{t - t>) + <^V)<(0^' - tj) (A.7) 
— Mr t(Mt) + b f Mt))(Mt') + b o Mt'))0(t ~ t') 



+ (Mt) + b Mt))(Mt') + b f Mt'))0(t' 
M£((b - b f ) [Mt)Mt')o(t - + MWitfW 
+ (6 - b f )b f <fo(t)<htf) + (/> bf (t)<p bf (t') 
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Imposing the boundary conditions tells us that 



Of 



. 



K^it') [(9, + «/)(rf W (t) + 6?^ W (t))*=vj = • ( A - 8 ) 

Hence if the linear combination of modes <j)j® = ^ + 6q r2 satisfies the boundary 
condition (d t + «o)^ (*)|t=*o = at i = i , i.e. 

Ko4 W (*o)+«t4 W (*0) 

and similarly for ij, the boundary conditions are obeyed. Finally is determined by the 
normalization condition (jL +a; 2 )G = — i<5. 

= (6 w _ 6 W) ^pww _ ' (A ' 10) 



J 2 ~ ^2 u, ri 

It is a standard exercise to show that the Wronskian (the Klein-Gordon inner product) is 
independent of t. 

The Lagrangian way treats all directions on the same footing. Recall that the frequency 
<jj is determined in terms of the eigenmodes of the spatial component of the Laplacian. 
The Lagrangian Green's function similarly uses eigenfunctions of the temporal Laplacian 
as building blocks, 

^ n \t) = -a 2 {n)^\t) . (A.ll) 



We assume that cr(0) = for simplicity. For each n there will be two real solutions to 

(n) (n) 

the eigenfunction equation <f>\ and 4> 2 ■ Together they form an orthogonal and complete 
set over a covering space containing the domain t £ [to,tf] (in general with the measure 



^m,n^i,j 



EKn)# ) (t)# ) (i') = ^4f ■ (A.12) 
i,n V-9(t) 

Here fi(n) is the measure on the 'dual' space. If the integral over t is over a non-compact 
domain, the sum over the eigenfunctions becomes an integral. Contrary to the statement 
below ( |A.3[ ), the boundary conditions for the Lagrangian Green's function are not satisfied 
by adding homogeneous terms. This can be traced back to the fact the delta function 
appearing in ( J^- + lv 2 )G = —id must obey the boundary conditions as well. This is done 
by the introduction of image charges in the covering space outside the domain t € [toify] 
of interest. An immediate consequence of the fact that the covering space is larger than 
the domain [ig, */] is that one expects the mode sum will be truncated to only those 
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modes which obey both boundary conditions. Here this is a direct consequence of the 
type of boundary conditions we are interested in. The conditions (dt + Ko,f)<j)\t=t ,t } = 
clearly leave the normalization undetermined. For a linear combination which satisfies one 
boundary condition, only a subset of the modes will obey the other. Choosing modes which 
manifestly obey the boundary condition at t = to, this subset is of course those modes for 
which 



W(t,) + W>dtW\t f ) = -k, (ff'{t f )+W>fi 

^ b (n) = «M n) (t/)+a^ W) (*/) _ b {n) (A 13) 

K f 4 n \tf) + a4 n \t,) f 

The Lagrangian Green's function then equals 

G%(*,o=iA^ e ( r l ■ (A - 14) 

The normalization A/g is determined by the condition (J™ + uj 2 )G = —id. For this one 
needs the explicit form of the eigenmodes. 

Because the boundary conditions uniquely determine the solution to a second order 
PDE, the Lagrangian Green's function eq. ( A.14j ) and Hamiltonian Green's function QA.7 ) 



are of course equal. In general this is difficult to show, but in specific cases one can do so. 
The most straightforward way is to decompose the Hamiltonian Green's function into the 
complete set of modes (f)& truno > _ if the domain x is non-compact and the mode sum in the 
Lagrangian Green's function becomes an integral, contour integration is another way to 
show equivalence. This contour will reveal multiple step functions, which contribute on the 
boundary, but are always constant on the domain [to, tt\. This way the Lagrangian Green's 
function recovers the statement that homogeneous terms enforce the boundary conditions. 

All these Green's functions can be written in a complex basis of eigenfunctions (p± = 
(pi ±i(p2 (Note the script <p notation for complex eigenfunctions). Sometimes this is a more 
condensed notation. One easily computes that 

«o<M*o) + d(j>i(t ) 



b bo =(t>i + bftfa , with &g 



^o<Mto) + d<h(t ) 



= + + • (A-15) 

The combination (1 ± ib^ ) can be related to the complex basis angle 6q : 

+ _ ko¥>+(*o) +dtp+(to) 



(1 - ibf) KQip- (t ) + dtp- (t ) 



. (A.16) 



Defining the complex analogue ip b = (p + + & c y?_, we find the relation between the real and 
complex bases, 

4>b = c KO (to)(p bo 
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= C KO (t )(p bo =C KO (t )b <p bo , 

c = («o¥>-(*o) + d<p-(t )) 

~ «W+Oo) + d(p + (t ) - (K cp-(t ) + d(f-(t )) ' 

Note that for real k the quantity 6g equals the inverse of its complex conjugate = 1/b . 
Hence Tp h = ip_ + b~ip + = \(tp+ + &</?-) = \<Pb- I n appendix |B| we shall argue that k 
should be treated as real throughout the calculation. For boundaries at a fixed time, i.e. 
initial conditions, k will be imaginary, but the correct results are only reproduced if one 
analytically continues from real k in the final correlation functions. We shall therefore 
treat k as real always. 

In this complex basis the Hamiltonian Green's function equals 



t 



KfH 

' ■ {ip bf (t)w t (ty(t-t>) + wMw f (t')e(t'-t)) . (A.i8) 



The normalization condition gives us that 

Mc = = — • (A.19) 

(1 - bibf) {<p+dip- - tp~dip+) 

The Lagrangian Green's function in the complex basis is 

(n) / ,x — (n) ,,n 

G L K ; 9 K0 (t, t) = iNt Kn f bo 2 ;, bo V ( A - 2 °) 



>)(t) + b^^\t)) U n \t>) + b^^\t>) 

-iA^JjM")- s/t 2 • 

Again to determine J\f£ one needs the explicit eigenfunctions: the terms with ^ (p(t)-tp(t')- : 
and l f{t)+ l P{t')+ generically do not vanish. If these terms solely contain contributions 
from the image charges, then A/jf ~ 1/(1 + bb). 

A.l Some simple examples 

In flat space R d_1,1 the spatial and temporal components of the kinetic operator sepa- 
rate cleanly and the complex mode functions are ip± = e ziztujt ; tp+ (t) = e mt . Therefore 
<p+d<p- - <f-d<P+ = -2iu), fdte i( - n - m ^ = 2ir5 n , m , and £ n ^ n ^~^ = 8{t - t'). Further- 
more 



b (n) = _ K ,f + in e2intoJ (A21) 



J o f 

kqj — in 
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Thus the Hamiltonian Green's function is 



G* a ™ (t, t') = e -^(*-*') - ^_^ e -M2t f -t-t') (A 22 



kq — iuj Kf+iujKQ — iuj J J 

The Lagrangian Green's function, on the other hand, gives 

trunc,K f in{t-t') _ Kg+in in(2tg-t-t') _ np-in in(t+t' -2t ) , J,n(t '-€) 
r<Lag i. ±l\ ■ ki-L 1 Kg— in n +in ' 

trunc,n f & i n (t-t?) _ Kg+in c in(2tn-t-t') 

= i(2N£)- lF4 • (A. 23) 

Recall the admonition below eq. ( [A. 17 ): k is assumed to be real. We see that for flat space 
= 1/2. The sum ranges over those modes for which 

2i(n 2 + KfK Q ) sin(2n(t/ - i )) - 2to(kj- - k q ) cos(2n(t/ - t )) = . (A.24) 

For kj = = Ki (hence boj = e nt °>-f), we indeed recognize the Green's function with 
Neumann boundary conditions at both ends: 

+ cos(o;(t + t' — to — tf) , 

trunc,Kr=0 ■ /, ,m . /„, t ,i\ 
1 ^.L J p in{t-t') I in(2t -t-t ') 

<3r~-=*5F ? £ — ^ — • (A - 25) 

where the sum is over the modes n = mir/2{tf — to), m £ Z. If we push the boundary 
tf off to infinity, the mode sum becomes an integral. Evaluating this integral by contour 
integration we recognize the step functions in the Hamiltonian Green's functions plus terms 
proportional to 9{2to — t — t') and 6(t + t' — 2to), each multiplying homogenous solutions of 
the kinetic operator. For the domain of interest t £ [to, oo), only the term with 9{t+t' — 2to) 
contributes, and we recover (part of) the homogeneous terms in the Hamiltonian Green's 
function. Choosing nt = oo = kq (hence boj = —e 2int °>f) we similarly recover the doubly 
Dirichlet Green's function. 

As we will show next (in appendix |B|) a particularly relevant choice of boundary cou- 
plings is K = —Kf = —itd (hence &o = bf = 0). From the expressions ( A.22j )-( |A.23j ), we 
see that this choice reproduces the flat Minkowski space Green's functions (subject to the 
mode selection rule ( A.24| ) reflecting the finiteness of the domain [io,i/])- Note that we 
obtained this result without an ie precription. This should be no surprise. The primary 
purpose of the ie prescription is precisely to ensure that the Green's function obeys the 
right boundary conditions. 
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B. Initial states in transition amplitudes, path integrals and fixed time- 
slice boundaries 



A naive Wick-rotation argues that the boundary action coupling constant k is imaginary 
for boundaries in time. For a real scalar field such a boundary condition, 

a^ = -*M0, (b.i) 

is at first sight inconsistent. We have stated that analyticity in the coupling constants 
for correlators computed in perturbation theory provides a resolution. One should treat 
k as real, and only analytically continue the correlation functions (including the Green's 
function) to complex or imaginary k. 

Here we show that this prescription is advocated by the relation of the path-integral 
to quantum-mechanical transition amplitudes. Recall that after a spatial Fourier trans- 
formation a field can be considered as an infinite set of harmonic oscillators, each with 
action 



S 



bulk 



r l f 


q 2 io 2 q 2 ' 


/ dt 




lt 


2 2 



(B.2) 



This action is obtained from the quantum-mechanical transition amplitude 



/ 



Vxe* sbM = (x N ,t f \e-^f-^\ Xl ,to) , H = ^+u 2 ^, (B.3) 



by splitting the interval tf — to into N smaller intervals of length (tt — to)/N, inserting 
N — 1 complete sets of \x) and N complete sets of \p) states, and taking the continuum 



limit N —* oo. This derivation makes clear that the action (B.2) has boundary conditions 
q(tf) = xn, q(to) = xi, and that the endpoints are not integrated over. Also clear is that 
temporal boundaries are quantum-mechanically on a very different footing than spatial 
boundaries. The latter simply affect the spatial modefunctions. Temporal boundaries, 
however, are encoded in the choice of initial and final state. 

For the free theory, a Gaussian integral, the exact answer for the transition amplitude is 
easily obtained. One substitutes the solution to the field equation with boundary conditions 
q(tf) = xn, q(to) = x\ into the action. Note that as the endpoints are not integrated over, 
the field equation is derived under the condition that the variation 5q vanishes on the 
boundary, Sq(tf) = 0, 5q(ti) = 0. One finds the well-known results (up to normalizations, 
which we ignore throughout this appendix) 



= De^ + c.c. , D- 



I 



2i sm(uj(tf — to)) 

Vx e° = exp 

_ JS ba < bulk {x N ,xi) 



( D 2 (e 2iujt f — e 2iujt °) 

-LO I C.C 



(B.4) 



Consider now the transition amplitude for a different initial state. In particular let us 
choose the harmonic oscillator vacuum |0) annihilated by d = ^ (ip + lux). This corresponds 
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to the Minkowski space vacuum for the field mode with frequency u. The transition 
amplitude {x N \e~ iH ^f- to ">\{)) can be obtained from the standard transition amplitude by 
the insertion of a complete set of states 

(xjv|e-^~* o) |0) = /dxi^le-^-^lxiXxilO) . (B.5) 



We can evaluate this expression in two ways. Either we can substitute the harmonic oscil- 
lator ground state wave function (xi|0) ~ e -ux\/2 anc j ^ e resu it (B.4) for the propagator. 
Performing the remaining Gaussian integral over x%, 



2 2 



dxi eiS »^ lHxN>xl)e ^ = g __iv (B 6) 

the result simply states that |0) is the zero energy eigenstate of the (normal-ordered) 
Hamiltonian. Or we can again derive a path-integral by splitting the interval tf — to into N 
smaller intervals, now inserting N complete sets of |x) and N complete sets of \p) states, 
and taking the continuum limit N — > oo. Doing so yields the bulk action ( |B.2j ) plus a 
boundary term 



Cjbulk+bdy 





' q 2 to 2 q 2 ' 


\ dt 


'to 


2 2 



^ ■ (B.7) 



As is clear from the ground state wave function (xi|0) the boundary coupling kq will 
be imaginary and equal to hq = —iuj. The Wick rotation intuition that the boundary 
couplings for spacelike boundaries are imaginary is confirmed. The answer for the transition 
amplitude (xjv"| e_i |0) ought then follow from solving the field equations for this 

action including the boundary term, and substituting the solution back. The extra insertion 
J dxi\xi) (xi\ means that the endpoint q(to) is now integrated over. The fluctuation 5q(to) 
therefore no longer vanishes and we obtain the field equations 

d 2 \ d 

W2 +u J 2 \ q(t) = , and - —q(t ) - K q(t ) = , (B.8) 

plus the implicit boundary condition q(tf) = xn- 

Because the coordinate q(t) is manifestly real, one has to give a prescription how to 



deal with the boundary condition (B.8) for imaginary k. It is quite obvious that insisting 
on q real, i.e. dq/dt(to) = = q(to), or insisting that the action remain real, q 2 — > \q\ 2 , will 
not reproduce the known answer (|B.6| ). However, if we simply proceed on the assumption 
that k is real, i.e. 



itut i v „—iijt\ 



q sol2 (t)=A(e wt + b Q e 
A(e^f + hoe-***) = x N , b = _^±^ e »«*o _ (B . 

Kq — liO 



the answer for the background value of the action, 



q bg,bulk+bdy _ luJ f X N A 2 h 2 -2iujt f \ ,p H| , 
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precisely reproduces the answer ( B.6| ) for kq = —iui (hence 6 = 0). This is therefore the 
prescription for dealing with imaginary boundary couplings: assume K is real until the final 
answer, and only then analytically continue. 

In the above example we have, of course, restricted ourselves to free field theory. 
One can repeat the whole exercise, however, with the inclusion of a bulk source term 
iS — > iS + J dtJ(t)q(t) representing interactions. Treating the source perturbatively, we 
expand into fluctuations £ around the background solution, q(t) = q so i(t)+^(t). Integrating 
the fluctuations out, we obtain for the action 



,-,2 



.2<f 



iJq 



+ Kf 



K - 



(B.ll) 



the result 



s b g ,m k+bdy{J . qsoi) = S b g;b uik +bdy{0) _ l j dtJ{t)qsol{t) 

- % - ! dtdt' J(t)G Kf;K0 (t,t')J(t') 



(B.12) 



where G KfiKo (t, t') is the Green's function of appendix [A]. Note that at endpoints where q(t) 
is not integrated over, i.e. when 5q{t en d) is constrained to vanish, £(i en d) also vanishes. At 
these points the Green's functions for the fluctuations £ therefore obeys Dirichlet boundary 
conditions with K en d = oo. For the transition function (x7v|e _J ^^ /_ * - ) |xi) we thus have 
Kf = oo = kq, whereas for the transition function (xAr|e _i ^^-'~ io )|0) we have Kf = oo, 
kq = —iu. Equivalence between the two transistion functions including bulk sources is 
thus established if 



/ 



dx\ exp 



iS^l k + M y{J-,q soh {x N , Xl ))\ (*i|0) 



exp 



■ cibq .bulk+bdy i T r w 

iS {J;qsoi 2 {XN)) 



K j — OO, 



.(B.13) 



The only dependence on x\ is in q so h{t) (eq. ( |B.4] )). Using the Green's functions of 
appendix which are derived with the assumption that k is analytic, it is an instructive 
exercise to verify that eq. ( B.13| ) is indeed true. The prescription to deal with imaginary 
k, by analytic continuation to imaginary values in the final correlation functions, therefore 
holds for perturbation theory as well. 

This example is an explicit manifestion of the fact that (in perturbation theory) all cor- 
relation functions are analytic in the coupling constants. This necessarily includes bound- 
ary couplings, which for a fixed time boundary correspond to initial conditions. 



C. Boundary field redefinitions in the presence of irrelevant operators 

We provide here the details behind the discontinuous shift of the field 4> on the bound- 
ary which effectively sets the coefficients of the relevant and irrelevant operators <f (f>d n <t>, 
§{d n (f)) 2 and § 4>d 2 (f) to zero. 
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After one integrates out high energy degrees of freedom, the most general form of the 
boundary action including the leading irrelevant boundary operators is 

Sbound = fd 3 X - - fyd n <i> - - ^dn4>d n 4> - . (C.l) 

Let us focus on the last operator for a moment. It is well known that in effective field 
theories (bulk) irrelevant operators of dimension p containing the factor df<p can be removed 
by a field redefinition at the expense of introducing irrelevant operators of dimension q > p 
p3| , |34|| , The only new element here is that the irrelevant operator is localized on the 
boundary. Generalizing, we see that the discontinuous field redefinition 

<f>(y)^Hy) + Hyo-y)jj<f>(y) (c.2) 

precisely generates a term that cancels the coefficient of j> (pd^cj) to first order in j3 c . To 
this same order the other couplings change as 

k' = k + Jji-™ 2 + 2k5 (°) + M<5'(0)) , 
^ = ^ + ^-2)5(0), 

0\\ = 0|| - Pc , 

(3' ± = (3 ± , (C.3) 



with primed quantities denoting the effective value after the field redefinition ( |C.2j ). Here 
m 2 is the bulk mass. We have ignored any bulk contributions to the boundary action of 
order <ft 3 and higher, having perturbation theory in mind. In appendix [S| we show that the 
explicit delta functions at zero argument, 5(0), serve to make all distributions conform to 
the boundary condition d n f(y) = —nf(y). 

To account for all couplings conflicting with the calculus of variations, fj,, /3m , and f3 c 
we combine the discontinuous field redefinition ( |C.2| ) with a discontinuous field redefinition 
of the form considered in section ^. 

0GO <Kv) + 8(yo - y) M(y) + u 2 d n <p( y ) + ...]+ 5( yo - y ) [a<p( y ) + ...] . (c.4) 

We have left the coefficient a arbitrary; as we will see there are additional compensations 
necessary beyond 5 = j3 c /M. Note that both a and a 2 have dimensions of M~ 1 . Consistent 
with the degree of approximation of the effective action, the field redefinition is an expansion 
in irrelevant terms to first order in M~ 1 . 

Formally we can solve for aj, 5 in terms of /i, P± and (3 C , so that the coefficients of 
the operators § (f)d n 4>, § (d n (j)) 2 and § 4>d 2 (j) vanish. The initial action Sbound of eq. flC.I|) is 
therefore equal to an effective boundary action 

a I j3 K e ff(ati,di) 2 

S eff = <b d d x '-tf (C.5) 



with the solutions for Qj substituted. (We absorbed the § (3ud l (j)di(f) into a momentum 
dependent n e ff(di)). At the end of the day we are only interested in the solution up to 
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linear order in /3±, (3 C . Higher order terms in f3± and (3 C would require the inclusion of 
higher order irrelevant operators for consistency. We may therefore linearize the problem 
and solve the system order by order in (3±, j3 c . Substituting the zeroth and first order terms 

ai = aw + a 12 — + 013 — + O(j^) , 



o 2 



a 



+ a 22 — + o 23 — + 0(^r) , 
+a 42 -+a 43 - + 0(^), 



(C.6) 



where aio is the solution given in eq. (|2.5| ), one finds (ignoring the zeroth order term in 

Sbound = S-i + Sq + S\ l + Si 2 , 



Su 



Si, 



So 



(fx cj)dl(j) 



M 



""10022 



«42 

2 



a 42 ai 



^«22(1 + — , 



Pc ( -«1QQ23 _ 043 

M\ 4 2 



0:43010 



A* n . «io- 

^«23(1 + — , 



1(1 + ^)2 

2 y 2 ' 



+ 



d xdnfidnC, 
Pc_ (_ a 23 

M \ 2 

d 3 x 05 n 



a 22 o 22 o 10 



M V 2 

"23010 A* 



A* /, . a io- 
-a 22 (l + — 



2 V 2 ; 



023(1 H 

4 v 2 



») 



M 



-012 



oi 2 aio . oi a 22 , « 
+ — n — 'H ) 



+042(1 + 



2 

3aio 



2 1 2 J 



)*(0) 



^«22 OlO, 

2 1 2 y 



-| ((ai 2 + 25 42 5(0))(1 + ^) - 022(1 + 010) 



+^ 
M 



-Ol3 

2 

K«23 



+aio(l + ^)*(0) 
013010 oi a 23 

n + n <K°) 



3flo 43 n . aio. 3oiq 

— 2 — (1 + — ) + «43(1 + — )5(0) 



2 (1 + ^) - f ((013 + 2a«*(0))(l + ^f) - 023(1 + aio) 



+ aio(l + ^(0) 



5-i 



/9_l / ai 2 oio 



M \ 2 

3i^04 2 aio 



xfn s . 42 2 . Oio, 

o(0) + — m (1 + — J + 



+- 



5(0) - « 



012 



9 + o 42 5(0))(l + } 



O42Ol0 ; - / , n s ~ x2( , n s 
5 (0) - O 4 2Oi Z (0) 

OiQ. 



(-(^ + 5 4 25(0))o 10 <5(0) + (-ai 2 5(0) + S 42 <5'(0))(1 + 



aio 
2 



-^ 2 (0) 



& O13O10 

Ml 2 



+— I — ~ — o(0) + — m (1 + — ) + 



aip, 
2 



043010 



y(0)-5 43 aio<5 2 (0) 
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+ g ^ ~ K VT~ O43<5(0))(l + — . 

-f (-(^ + S 4 35(0))a 10 5(0) + (-a 13 <5(0) + 5 43 <J / (0))(1 + ^) 
+^(1 + ^K(0) 



(C.7) 



These equations can be explicitly solved (e.g. 023 = 0). For the case \i = (as in Bunch- 
Davies for instance), and hence a\o = 0, the solutions are easily found: 023 = 0, 022 = 
— 1, 5?42 = 0, S43 = —1, «i3 = (25(0) — 3H), a\2 = K, with the answer for S-i: 



d A x 



K+ M K 



M 



(m 2 + k 2 - 3Hk) + 4rc|p(0) 



(C.8) 



As we will show in the next appendix, the term k(3 c 5(0) solely served to make all distri- 
butions consistent with the boundary condition d n f(y) = —nf(y). We may therefore drop 
this term, as long as we remember this. 



D. Distributions, boundary conditions and the equivalence between per- 
turbation field theory and field redefinitions 



The result ( C.8 ) for the effective boundary action after the field redefinition suggests that 
the correction to the two-point function due to irrelevant operators contains delta func- 
tions at zero argument. In the perturbative Feynman diagram approach to the two-point 
function, which we perform in the next appendix, we shall find no explicit 5(0) terms. Yet 
the two approaches are manifestly equivalent, so somehow a further step is needed in the 
field redefinition approach to explain why no 5(0) term arises in the two-point correlator. 

To understand the equivalence between the two approaches better, consider a matrix 
integral simplification of the path integral 

A x-j k, x^ x^ (3 - ■ x* 1 x^ 



(x k x l ) = J dx i x k x L e "2 "2 "a . (D.l) 

Here Hij and /3ij correspond to the boundary interactions; whereas A{j is the kinetic oper- 
ator. We will now evaluate this integral in two ways (1) by a saddlepoint approximation, 
i.e. a Feynman diagram expansion with 0ij treated as an interaction, and (2) by a field 
redefinition which absorbs at the expense of redefining mj. Expanding the answer (2) 
to linear order in we should reobtain the Feynman diagram result. 

The Feynman diagram approach: Expanding to linear order in /3y we find that 
<*V) = I dxW ( 1 - AA + ..) e" 1 ^ 2 



2 



dJ k dJi V 2 dJi dJj 



(D.2) 

j=o 
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Here N is an unimportant normalization, and we have introduced the Green's function 
G K = (A + /i)~ 127 . The Gaussian integrals are easily evaluated to 



(x k x l ) = G kl - Pij&jGf + (1 - ^Gjf . (D.3) 
We clearly recognize the connected and loop diagrams. 

Field redefinitions: The field redefinition is designed such that to first order the contri- 
bution from the kinetic part cancels the fiij factor. Hence 

X «_ I <_^* X *= (i-9l\ x = Sx , (D.4) 

where we have introduced a second Green's function G = A -1 . Note that G ^ G K . Under 
this field redefinition the integral becomes 

{x k x l ) = SpSg / dx\3&c\x p x q e 2 . (D.5) 

The Jacobian |Jac| will contain the loop diagrams. Our interest only extends to connected 
diagrams and we may therefore ignore it. Expanding to linear order in /3y we get 

(x k x l ) = S k S l q J dx\3ac\xPx q e~^ T ^ A+K ~^ L ~ !L ^~) X . (D.6) 

The term proportional to k/3 is exactly the problematic one, as we will see. Thus the 
two-point function is easily evaluated to 



{x k x l ) = S k p S\ [G K + G K ( ^Y^ ) G k ) 



pq 



Gk-^y [{A + k)G0 - kGP - P t Gk + (3 T G(A + k) ) G H 



ki 



= (G K - G R f3G K ) kl . (D.8) 

In the second to last step we recalled that G K = (A + k)^ 1 . We see that we exactly 
reproduce the connected diagrams as expected. 

Applying these lessons to field redefinitions on the boundary: As eq. ( p.7| ) shows 
field redefinitions which are localized on the boundary ought to have no effect on bulk 
correlators. This means that the fourth term in (D.7) ought to reproduce the Feynman 
diagram computation. If the kG/3 factor contains the (5(0) term, this appears not to be the 
case. The resolution follows from repeating the steps (p.7|) in detail in the field theory. 



27 In the following we will use matrix and index notation interchangeably. It should be clear from the 
context which notation is being used. 
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If we consider the index % as the location in the y direction, we easily see that in 
the (free) field theory of section || with boundary interaction ( |C.1| ) and (3± = /3m = the 
matrices A, k, (3 correspond to the differential operators. 

Aij = A(y 1 ,y 2 ) = <%i - y 2 )did 2 , 

k = K co 5(y - yi)<%i - V2) , 

P = PcoKvi - V2)Kyo - yi)n 2 . (d.9) 

We have given the couplings a subscript co to distinguish them from the matrix operators. 
We easily compute that 

A- 1 (y 1 ,y 2 ) = G(y u y 2 ) with UG = -<%i - y 2 ) and d y G = 
G(3 = /3 co J dy 2 G( yi ,y 2 )D 2 5(y 2 - ys)5(yo -2/2) 

= -Mm - ys)%o - ya) ■ (d.io) 

The transformation £ is thus indeed the one we consider in eq. (C.2). 

dy 2 S( yi ,y 2 )<p(y 2 ) = / dy 2 S(y 1 - y 2 )<j>( y2 ) + -^(yi - ^2)^(^0 - V2)(t>{y2) 

= Hyi) + ^(y -y 1 )<P(y ) . (D.ll) 

We also see that 

= -K co /3 co 5(yo - yi)%i - ^2)^(^0 - 2/2) ( D -12) 
and hence that it contains the problematic 5(0) term: 

dyidy 2 4>{ yi )4>(y 2 ) [kG0\ (y la y 2 ) = -K co /? co ^(y ) 2 5(0) . (D.13) 



Because this 5(0) term is present in the action, we expect it to be present in the two-point 
correlator as well. Indeed a straightforward computation gives (Note that A' 1 = G li : the 
Green's function obeying the boundary condition d y G = —kG.) 

(<t>( yi )<j>(y 2 )) = G K ( V1 , y 2 ) + 5(y x - y ) Gk{iJ °> m) Pco + fl, Gr{v ^ Vo) 6(y - y 2 ) 

-G K ( yi ,y )6(0)G K (y ,y 2 ) . (D.14) 

The Feynman diagram computation, however, has no 5(0). How can this agree? We have 
seen the explicit steps we need to do to get the Feynman diagram answer. Surprisingly 
when we implement them here, the 5(0) cancels. We will need that 

A + k = 5(y± - y 2 )did 2 + K co 5(y - yi)5(y 1 - y 2 ) 

- ~5(yi - y 2 )Oi + 5(yi - y 2 )5(y - yi)(di - k) . (D.15) 

Its inverse, the Green's function G K , obeys a slightly different differential equation, however. 
As we also discussed in appendix [A], acting with the Laplacian on G K returns the delta 
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function 5 K (y\ — 2/2) in the space of functions obeying d y f(yo) = —Kf(yo). There are 
additional contributions from image charges which guarantee that on the boundary dy5 K — 
— kS k . Now repeating the steps from eq. ( p. 7 ) 



(0(2/1)0(2/2)) = (2/1,2/2) 

+ / dy 3 dy4(3coG K (yi, y 3 ) (-D 3 + K co 5(y - 2/3) )%3 ~ 2/4)<%4 - 2/o)^^r-^A: 



+ J dy 3 dy A (3 co GK ^ V ^ 5(y - ^ 4 )(— D 3 + K co 5(y - y 3 ))5(y A - y 3 )G K (y 3 ,y 2 ) 
- G K (2/1 , 2/0 ) 5 (0) G K (y , y 2 ) 

= G K (y 1 ,y 2 ) - ^-5 K (y x - yo)G K (y ,y 2 ) - ^-5 K (y 2 - 2/0)^(2/1,2/0) , (D.16) 

we recognize that the sole function of the 5(0) term in the action is to correctly implement 
the boundary conditions for the transformation S = 1 — G/3/2. Indeed it is clear from 
the matrix analogue that had we started with a transformation S K = (1 — G K /3/2) no 
distributions at zero argument 5(0) would have been generated at all. 

The distribution 5 K (yi — 2/0) with the correct boundary conditions which thus appears, 
has no support deep in the bulk, 5 K (y± — 2/0) vanishes for 2/1 3> 2/0 of course. The lesson we 
extract from this exercise is that for bulk correlators we may ignore the 5(0) term in the 
action. 

Other field redefinitions: If field redefinitions <p(y) — * <p(y) + 5(y — yo)a(j)(yo) to remove 
the irrelevant operator <f> fl co (\)d\§ leave no trace in bulk correlators, an obvious question is 
why the "theta" transformations do contribute. They do, and why follows from repeating 
the above steps for that case. Consider for simplicity only the relevant correction fj,. In 
the above language it corresponds to choosing 

Pn = HcoKvi ~ 2/o)%i - y2)d yi ■ (D.17) 

Therefore 

Gf3^ = j dy 2 d y2 G(y 1 ,y 2 )5(y 2 - y )5(y 2 - y 3 ) = d V3 G(yi,y 3 )5(y 3 - y ) . (D.18) 

Now (in Minkowski space) we can show that this is exactly the step function transformation. 
Upon use of the identity 

G/3 M = -d yi G(y x ,y 2 )5(y 2 - 2/0) , (D.19) 

we can take one more derivative to obtain 

-dlG( yi ,y 2 )5(y 2 - y ) = -5( Vl - y 2 )5(y 2 - y ) . (D.20) 

Thus Gfin precisely has the property that it's derivative is the delta function — it is 
therefore proportional to the 9 function. G, moreover, is the Neumann Green's function. 
Hence Gj3 is zero in the bulk, it is precisely equal to 9(yo — y\). 
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Because 9(i/q — y{) is of measure zero, the statement that the second and third terms 
G(5^G~ l arising from the explicit field redefinition do not contribute to the bulk, is now 
manifest. Thus the fourth term — the one that comes directly from the action — ought 
to reproduce the Feynman diagram result. Indeed it is easy to see that 

J dy 2 dy 5 G K (y 1 , y 2 )(-K co Li co 5(y2 - yo)dG(y 2 , ys)<%5 - ya))G K (y 5 , y 4 ) (D.21) 

precisely reproduces the perturbative Feynman diagram calculation, when we use the just 
derived result that dG(y ,y ) = 0(0) = 1/2. 



E. Power spectrum corrections from perturbation theory 

Aside from using field redefinitions on the boundary, one can also use field theory pertur- 
bation theory to compute the corrections to the power spectrum. For completeness we give 
that calculation here. The answer is, of course, the same as in eq. ( 4,19| ) to first order in 
Pi- 

The first order correction to the (connected) two-point correlation function by a gen- 
eralized two-point vertex 

= - [ d A x^-gh 2 = - [ d 4 x 3 d 4 x A 73^) y/r^) *(S3, *4) 0(a . 3)0(a . 4) (E.!) 



2 



IS 



K/ (0(xi)0(x 2 )) K = -i J d 4 x 3 d 4 X4y / -g(x 3 )v / -g(x 4: ) X(x 3 , x i )G K (x 1 ,x 3 )G K (x 2 , x 4 ) . (E.2) 

Here ft/, k denote the future-out and past-in state and G K (x\,x 2 ) is therefore the Green's 
function satisfying (Di — m 2 )G K (x\, x 2 ) = i5 4 (x\ — x 2 ) / —g with the boundary conditions 

Og d T j 1 G K (xi, X 2 )\fj 1=r j = — kG k (x\, %2)\r]i=r)o 

lim a~ l d m G K (xi, x 2 )\ m=r]f = -KtG R (xi, x 2 )\ Vl=ri . (E.3) 

For the boundary interaction due to the leading boundary irrelevant operators ( |4.6| ), the 
(spacetime dependent) coupling X(x 3 ,X4,) considered as a derivative operator equals 

A(* 3 ,, 4 ) = 2 ( ^3-^0)^(x3-X 4 )^3 -„ 4 ) \ x 

V a{r] 3 ) 0^(7/3) 0(773) ) 
+ zJkl7 d vs d m + ttJ^tj (D m d m + D m d m ) + (d m + d m ) 



a?{r, 3 )M 1 ' a 2 (r] 3 )M m m ' 2a 2 (%)M v "* m ' m mJ ' 2a(r] 3 ) 

(We purposely avoid integrating by parts, because A(#3,x 4 ) arises from a boundary action 
rewritten as bulk interactions. Integration by parts would make this origin less clear.) 
Inserting this expression and the appropriate FRW quantities in eq. ( [E.2D , we obtain after 
a spatial Fourier transform 
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2i d> (Fx^d? xi(fx 2 age 

J tl3=V4=V0 

2M 



3H 



^3 ^4 



O 1 " "0 

/i 

2a 

2i(2vr) 3 5 3 (A?i + fc 2 )ao x 



+—(^3+^4 



/3 C 



alM + a^M^ 4 + M 



k 2 

□ 3 - -D 4 - 7r-(d m + d m ) \ 



3H i 

2a 



G K (h,r]i,r) 3 )G K (ki,r]2,ri4) 



V3=V4=110 



-2i(2ir) 3 5 3 (ki + fe 2 )ao x 



2a 4 



fc 2 

m 2 - -| + 3ffK 



a, 



o 



/UK 



G K (k 1 ,rj 1 ,r] )G K (k 1 ,r]2,r]o) 



(E.5) 



In the first step we related the double normal derivative D n d n to the Laplacian. In the 
second step we used both defining property of the Green's function (□ — m 2 )G K = i5f./a 4 
and the boundary condition d v G K = —ciokG k . Recall the expression for G K (k\, 771,772) from 
eq. (fT?]) in terms of the basis functions <fidS,±- For the power spectrum we are interested 
in the equal time two-point correlator for 771 = 772 — > 0. In that limit the Green's function 
only retains the retarded contribution 



(E.6) 



Below we shall see that for the inflationary power spectrum, we should choose Kf = k. The 
equal-time correlator at 771 = 772 — > therefore equals 



lim K (</)(rii,h)<t>(vl,h)) K = -2i(2vr) 3 5 3 (A: 1 + k 2 )a&£ K (mM K (Vo) x 

J71-+0 

~k 2 {P\\- Pc) , k 2 P ± p c m 2 3(3 C H. 

+ — —-K(fx-——) 



M 



M 



M 



(E.7) 



Using the proportionality relation ( f4,13| ) between the basis-functions ipb and ip b as 77 — > 0, 
plus the expression for the zeroth order two-point correlator we obtain 



lim . «;<0(r?i, A;i)0(r?i, fc 2 )>« = (4> )o 

?71— »0 



1 



b-1 

+ 



-2ial)^ K (r] ) 2 x 



k 2 (P\i - 0c) k 2 (3 ± (3 c m 2 



M 



M 



M ' 



(Ei 
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Finally substituting the explicit expressions for (pb K , 



lim . K (4>(r] 1 ,k 1 )4>(r] 1 ,k 2 ))K = {<t> )o 

?7l— »0 



1 



-2i7T— 2 



6 - 1 / 4iJ 

/9c) k 2 /3 ± /3 c m 2 



+ 



M 



M 



kr] )x 



M ' 



. (E.9) 



with the obvious shorthand H^ v = H u + bH u . 

The power spectrum of inflationary density perturbations due to spontaneous pair 
production in a gravitational background is obtained by the optical theorem from the 
two-particle cut of the one-loop vacuum amplitude (k\k). 



Pd 6 k 



(47T) 



lim Im 



i (<f>{m,h)<t>( r ll,k2))K \ d\k\ 



(27r) 3 J/1-+0 



(E.IO) 



This shows that Kt = R. (Note the factor of i; this is a consequence of our normalization 
for the Green's function.) The imaginary part of the (Feynman time-ordered) Green's 
function is also known as the Wightman function. In contrast to the Green's function, the 
latter is a homogeneous solution to the field equation. We thus find 



P. 



K. + S0 



P.. 



(E.11) 
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1 - b\ H b J-krjo) 
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which agrees with eq. ( [4,19| ). 
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